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ABSTRACT 


Ocean  wave  eharacteristies  in  the  western  Atlantic  Ocean  (Hurricane  Region)  to 
tropical  cyclones  have  been  well  identified,  but  not  the  regional  seas  in  the  western 
Paeifie  sueh  as  the  South  China  Sea  (Typhoon  Region).  This  is  due  to  the  laek  of 
observational  and  modeling  studies  in  the  regional  seas  of  the  western  Paeifie.  To  fill  this 
gap,  Wavewatch-III  (WW3)  is  used  to  study  the  response  of  the  South  China  Sea  (SCS) 
to  Typhoon  Muifa  (2004).  The  major  purposes  are  to  find  the  similarity  and  dissimilarity 
of  wave  characteristics  between  the  two  regions,  and  to  evaluate  the  WW3  capability  to 
typhoon  forcing.  The  WW3  model  is  integrated  from  the  JONSWAP  wave  speetra  with  a 
tropieal  cyelone  wind  profile  model,  simulating  Typhoon  Muifa,  from  OOOOUTC  16 
November  to  1200UTC  25  November  2004.  Since  TY  Muifa  entered  the  SCS  as  late  as 
19  November,  the  model  eomputation  of  the  first  three  days,  from  16  November  to  18 
November,  could  be  eonsidered  as  the  ‘spin  up’  period  of  WW3  model.  This  study 
shows  strong  similarities  in  the  responses  between  Hurrieane  and  Typhoon  Regions, 
ineluding  strong  asymmetry  in  the  signifieant  wave  height  (H^)  along  the  typhoon 

translation  traek  with  the  maximum  in  the  right-front  quadrant  of  the  typhoon  center, 

and  asymmetry  in  the  direetional  wave  speetra  at  different  loeations  (frontward, 
baekward,  rightward,  and  leftward)  around  the  typhoon  eenter.  The  unique  features  of 
the  SCS  wave  eharaeteristies  to  Muifa  are  also  discussed. 
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I.  INTRODUCTION 


A  moving  tropical  cyclone  is  an  intense  source  of  surface  wind  stress  that  causes 
many  signifieant  ehanges  in  oeean  wave  eharaeteristies  sueh  as  signifieant  wave  height, 
directional  wave  spectra,  and  wave  propagation.  These  features  have  been  well  identified 
in  open  oeeans  and  the  western  Atlantie/eastern  Paeific  regional  seas,  that  is,  the 
Hurrieane  Region  (hereafter  referred  to  HR).  The  tropieal  eyelone  foreed  eomplex  wave 
field  is  very  important  to  the  wind-wave  interaetion.  A  hurrieane  with  intense  and  fast- 
varying  winds  produees  a  severe  and  eomplex  oeean  wave  field  that  ean  propagate  for 
thousands  of  kilometers  away  from  the  storm  eenter,  resulting  in  dramatie  variation  of  the 
wave  field  in  spaee  and  time  (Barber  and  Ursell  1948).  To  investigate  the  wave 
eharaeteristies,  the  direetional  speetra  of  hurrieane  generated  waves  were  measured  using 
various  instruments.  For  example,  the  feteh  effect  was  detected  in  the  Celtie  Sea  using 
the  high-frequeney  radar.  The  wave  eharaeteristies  were  obtained  for  the  northeastern 
Paeifie  during  passage  of  storm  using  the  synthetie  aperture  radar  image  from  ERS-1 
satellite  (Holt  et  al.  1998).  The  spatial  wave  variation  of  hurrieane  direetional  wave 
speetra  was  identified  for  both  open  oeean  and  landfall  eases  using  the  US  National 
Aeronautics  and  Space  Administration  (NASA)  seanning  radar  altimeter  (Wright  et  al. 
2001;  Walsh  et  al.  2002). 

The  oeean  wave  response  identified  in  HR  is  a  signifieant  right-forward-quadrant 
bias  in  the  signifieant  wave  height.  During  the  passage  of  Hurrieane  Bonnie  (1998)  in  the 
Atlantie  Oeean,  both  observational  (Wright  et  al.  2001)  and  modeling  (Moon  et  al.  2003) 
studies  show  that  the  significant  wave  height  reaehes  14  m  in  the  open  oeean.  The 
maximum  wave  heights  appear  in  the  right  forward  quadrant  of  the  hurricane  eenter  and 
propagate  in  the  same  direetion  as  the  hurrieane.  Moon  et  al.  (2003)  simulated  the  wave 
eharaeteristies  sueeessfully  using  the  wave  model  WAVEWATCH-III  (hereafter  referred 
to  WW3)  and  found  that  the  hurrieane-generated  wave  field  is  mostly  determined  by  two 
faetors:  the  distanee  from  the  hurrieane  eenter  or  radius  of  maximum  wind  and  hurrieane 
translation  speed.  Eor  the  ease  of  a  hurrieane  with  low  translation  speed,  the  dominant 
wave  direetion  is  mainly  determined  by  the  distanee  from  the  hurrieane  eenter. 
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Most  of  the  observational  and  modeling  studies  on  oeean  waves  generated  by 
tropical  cyclones  are  concentrated  in  HR.  Few  observational  and/or  modeling  studies 
have  been  done  in  the  Typhoon  Region  (hereafter  referred  to  TR),  especially  in  the  South 
China  Sea  (SCS).  The  SCS  is  one  of  the  largest  marginal  seas  of  the  Western  Pacific 
Ocean,  extending  across  both  tropical  and  subtropical  zones  and  encompasses  a  total 
surface  area  of  3.5xlO*’A:m^  (Figure  1).  Also  due  to  its  semi-enclosed  nature,  the  SCS  is 
subject  to  high  spatial  and  temporal  variability  from  external  forcing  factors.  One 
significant  source  of  the  SCS  variability  is  the  tropical  cyclones  that  routinely  affect  the 
region.  The  WW3  has  been  implemented  and  verified  for  the  SCS  using  the 
TOPEX/Poseidon  (T/P)  satellite  data  (Chu  et  al.  2004).  However,  there  is  no  modeling 
study  on  SCS  waves  to  typhoon  winds.  Our  goal  in  this  research  is  to  identify  if  those 
effects  occurring  in  HR  still  exist  in  TR.  More  specifically,  we  study  the  SCS  responses 
to  Typhoon  Muifa  (2004)  using  WW3,  which  was  forced  by  a  high-resolution  wind  field 
computed  by  a  Tropical  Cyclone  Wind  Profile  Model  (TCWPM)  proposed  by  Carr  and 
Elsberry  (1997). 

The  outline  of  this  thesis  is  as  follows.  The  geography  and  climatology  of  the 
SCS  are  described  in  Chapter  II.  The  descriptions  of  the  SCS  tropical  cyclones  and 
Typhoon  Muifa  (2004)  are  in  Chapter  III.  Chapter  IV  discusses  the  processing  of  high 
resolution  typhoon  winds.  The  model  features  of  WW3  and  its  implementation  are 
described  in  Chapter  V.  This  study  also  reviews  several  important  recent  researches  in 
Chapter  VI.  The  numerical  simulation  and  the  wave  characteristic  of  the  SCS  are 
analyzed  and  discussed  in  Chapter  VIE  The  conclusion  and  some  suggestions  are 
provided  in  Chapter  VIII. 


2 


China 


Longitude 


Figure  1  Geography  of  the  South  China  Sea  (from  Chu  et  al.  2004). 
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II.  SOUTH  CHINA  SEA 


The  SCS  is  one  of  the  most  important  marginal  seas  of  the  Western  North  Paeifie 
Oeean.  It  is  subjeeted  to  seasonal  monsoon  foreing  with  the  northeasterly  in  the  winter 
and  the  southwesterly  in  the  summer.  Tropieal  eyelones  also  invade  the  SCS  in  any 
month  of  a  year. 

A,  GEOGRAPHY 

The  SCS  is  loeated  among  the  Asian  landmass  to  the  west  and  to  the  north, 
Taiwan  to  the  north-east,  the  Philippine  Islands  to  the  east,  and  Borneo  to  the  south 
(Figure  1).  The  SCS  region  extends  from  the  equator  to  23°N  and  from  99°E  to  121  °E 
and  the  total  surfaee  area  is  about3.5xlO®A:m^.  The  SCS  eonneets  to  the  East  China  Sea 
through  the  Taiwan  Strait,  to  the  Paeifie  Oeean  through  the  Luzon  Strait,  to  the  Sulu  Sea 
through  the  Balabae  Strait,  to  the  Java  Sea  through  the  Karimata  Strait,  and  to  the  Indian 
Oeean  through  the  Strait  of  Malaeea.  All  of  these  straits  are  shallow  and  narrow  exeept 
the  Luzon  Strait  with  the  maximum  depth  about  1800  m. 

The  SCS  is  a  semi-enelosed  tropieal  sea  with  a  large  abyssal  basin,  of  whieh  the 
maximum  depth  reaehing  as  far  as  5000  m  at  the  geographie  eenter  (Su  1998).  The 
ellipse-shaped  basin  is  about  1900  km  along  its  major  axis  (northeast-southwest)  and 
approximately  1100  km  along  its  minor  axis.  Wide  eontinental  shelves  appear  in  the 
northwest  and  southwest  of  the  basin  and  steep  slopes  in  the  eentral  portion.  Many  reef 
islands  and  underwater  plateaus  are  seattered  throughout  the  SCS.  Lrom  the  Taiwan 
Strait  to  the  Gulf  of  Tonkin,  the  eontinental  shelves  are  about  70  m  deep  and  150  km 
wide.  In  the  south  end,  the  Sundra  Shelf  is  the  submerged  eonneetion  between  the 
Southeast  Asia,  Malaysia,  Sumatra,  Java,  and  Borneo,  and  is  100  m  in  depth  around  its 
eenter  (Li  1994). 

B,  CLIMATOLOGY 

In  the  winter  season,  winds  are  generally  out  of  the  northeast  and  relatively 
strong,  while  in  the  summer  season,  the  winds  completely  reverse  and  become 
southwesterly  and  weak.  Lurthermore,  synoptic  systems  often  pass  by  the  SCS  and  cause 
temporally  and  spatially  varying  wind  field.  Tropical  cyclones  can  pass  through  and  even 

generate  within  the  SCS  in  any  month  of  the  year. 
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1.  Monsoon  Winds 

The  seasonal  movement  of  the  equatorial  pressure  trough,  also  known  as  the 
intertropieal  eonvergenee  zone  (ITCZ),  produees  seasonal  wind  flow  over  the  SCS  (Gray 
1968).  A  elassic  definition  of  the  monsoon  ean  be  referred  to  Ramage  (1971).  The 
monsoon  area  is  an  encompassing  region  with  January  or  July  surface  circulation,  in 
which  the  prevailing  wind  direction  shifts  by  at  least  120°  between  January  and  July,  the 
average  frequency  of  prevailing  wind  directions  in  January  and  July  exceeds  40  %,  and 
the  mean  resultant  winds  in  at  least  on  of  the  months  exceed  3  m/s. 

The  SCS  satisfies  all  above  features  and  experiences  both  winter  and  summer 
monsoons  every  year.  The  winter  monsoon  is  known  as  the  ‘northeast  monsoon’  and  the 
summer  monsoon  as  the  ‘southwest  monsoon’.  These  names  are  derived  from  the  low- 
level  prevailing  winds  of  the  two  seasons  (Figure  2).  The  winter  monsoon  onsets  in 
around  November  and  retreats  in  late  March  while  the  summer  monsoon  onsets  in  late 
May  and  retreats  in  late  September.  Between  them  are  the  brief  transitional  periods 
(Cheang  1987). 

2.  Summer  Monsoon 

In  July  and  August,  temperatures  from  the  Asian  landmass  to  the  north  of  the  SCS 
reach  the  annual  maximums  and  produce  a  lower  pressure  over  the  continental  region.  In 
contrast,  cooler  air  over  the  SCS  and  the  southward  oceanic  region  produces  a  higher 
pressure.  An  equatorial  trough  over  the  central  Philippines  extends  northwestward  to  the 
low  pressure  over  the  Tibetan  Plateau.  The  pressure  gradient  between  the  warm 
continental  regions  (north)  and  the  cooler  oceanic  regions  (south)  causes  air  to  flow 
northwestward  south  of  the  equator.  The  air  flow  then  turns  to  flow  northeastward  as  it 
crosses  the  equator  and  becomes  the  summer  monsoon.  The  pressure  gradient  during  the 
summer  monsoon  season  is  relatively  weak,  and  produces  southwest  monsoon  winds 
(about  3  m/s)  over  the  SCS  (Ramage  1971). 

3.  Winter  Monsoon 

The  summer  southwest  monsoon  begins  to  retreat  in  September,  as  the  Asian 
landmass  begins  to  cool  down  and  a  high  pressure  starts  to  build  over  this  region.  Air 
temperatures  over  the  oceanic  region  remain  warm  and  the  pressure  gradient  begins  to 
reverse.  In  October,  the  equatorial  trough  begins  to  move  rapidly  southward.  By  the 
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mid-October,  the  trough  is  along  a  line  from  the  center  of  the  Bay  of  Bengal  to  the  north 
eoast  of  New  Guinea.  North  of  the  trough,  northerly  winds  prevail  as  a  high  pressure 
centered  over  the  eentral  Asian  landmass  and  eontinue  to  build,  while  in  the  south  side  of 
the  trough,  the  southeast  monsoon  still  dominates  (Ramage  1971). 

The  winter  northeast  monsoon  begins  to  set  up  in  November  when  the  equatorial 
trough  moves  south  of  the  equator.  The  monsoon  intensifies  over  the  SCS,  and  monthly 
mean  wind  speed  inereases  to  7  m/s.  By  December,  the  equatorial  trough  is  loeated  north 
of  Australia  (5°  S).  The  high  pressure  is  firmly  established  over  the  Asian  landmass  and 
intensifies  the  pressure  gradient  between  the  continental  and  the  oeeanic  regions.  The 
northeast  monsoon  reaehes  its  strongest  intensity  over  the  SCS  (8  to  10  m/s).  The  flow 
over  the  SCS  is  northerly  to  northeasterly  winds  north  of  the  equator.  South  of  the 
equator,  the  reversal  in  the  sign  of  the  Coriolis  foree  eauses  the  flow  to  turn  eastward  and 
becomes  northwesterly  to  westerly  (Ramage  1971). 

The  northeast  monsoon  continues  over  the  region  until  April  when  the 
temperatures  over  the  Asian  eontinent  start  to  inerease,  and  when  the  equatorial  trough 
begins  to  move  to  the  north.  Winds  in  the  northern  SCS  remain  northeasterly,  but  are 
weakened.  In  May,  the  northeast  monsoon  eompletely  disappears,  as  the  Asian  eontinent 
continues  to  warm  and  the  pressure  gradient  between  the  continental  and  oceanic  regions 
reverses  (Ramage  1971). 
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Figure  2  The  SCS  monsoon  wind  pattern,  (a)  Summer  southwest  monsoon,  and  (b)  winter 

northeast  monsoon  (from  Cheang  1987). 
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III.  TROPICAL  CYCLONES 


A.  BACKGROUND 

Figure  3  shows  the  spatial  distribution  of  total  number  of  tropieal  eyelones 
passing  by  from  1945  to  1988  (Neumann  1993).  During  that  period,  about  26  tropical 
cyclones  are  generated  every  year  over  the  Western  North  Pacific  Ocean.  Among  them, 
quite  a  few  pass  by  the  SCS.  Besides,  a  tropical  cyclone  can  pass  by  the  SCS  in  all  the 
months  within  a  year  (McGride  1995).  The  primary  reason  is  the  persistently  warm  sea 
surface  temperature  and  the  location  of  the  intertropical  convergence  zone  (ITCZ). 
Annual  reports  on  tropical  cyclone  by  the  United  States  Navy  Joint  Typhoon  Warning 
Center  (JTWC)  show  that  about  80  %  of  the  tropical  cyclones  are  initially  generated  in 
the  monsoon  trough  (McGride  1995). 

Figure  4  shows  the  frequency  of  tropical  cyclone  occurrence  in  the  Western  North 
Pacific  Ocean  from  1945  to  1988.  The  solid  line  indicates  the  tropical  cyclones  with  the 
maximum  wind  speed  higher  than  17  m/s,  while  the  shadow  area  is  for  the  maximum 
wind  speed  higher  than  33  m/s. 

B,  TYPHOON  MUIFA  (2004) 

Typhoon  (TY)  Muifa  is  one  of  the  four  tropical  cyclones  passing  by  the  SCS  in 
2004.  It  was  formed  on  1 1  November  and  weakened  over  land  on  26  November.  The 
best  track  passage  of  TY  Muifa  (Figure  5)  and  best  track  record  (Table  1)  were  provided 
by  the  JTWC  (2005). 

1.  Forming  over  the  Western  Pacific 

TY  Muifa  was  first  originated  as  a  tropical  depression  on  1 1  November  2004  in 
the  Western  North  Pacific  Ocean.  It  moved  steadily  northwestward  passing  north  of 
Palau  before  entering  the  Philippine  Sea.  The  disturbance  was  first  mentioned  by  JTWC 
at  1600UTC  13  November,  approximately  400  km  north  of  Palau.  The  depression  was 
developed  into  a  tropical  storm  on  14  November.  According  to  its  strength,  the  Japan 
Meteorological  Agency  (JMA)  first  named  it  Muifa.  However,  the  Philippine 
Atmospheric,  Geophysical,  and  Astronomical  Services  Administration  (PAGASA) 
assigned  the  name  Unding  at  OOOOUTC  14  November  after  the  tropical  cyclone  had 
invaded  their  area  of  responsibility. 
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Tropical  Storm  (TS)  Muifa  was  steadily  upgraded  and  moving  west- 
northwestward  to  the  Philippines.  It  began  the  eloekwise  loop  at  1200UTC  17  November 
and  eontinued  for  several  days.  On  18  November,  the  intensity  of  TY  Muifa  reaehed  the 
maximum  intensity  at  54.0  m/s  (115  knots)  at  1200UTC  and  then  began  a  remarkable 
weakening  phase. 

The  satellite  image  of  TY  Muifa  (Figure  6)  was  obtained  from  the  Tropieal 
Rainfall  Measuring  Mission  (TRMM)  satellite.  The  TRMM  is  a  joint  mission  between 
NASA  and  the  Japan  Aerospaee  Exploration  Ageney  (JAXA)  and  eapable  with  both 
passive  and  aetive  sensor  array.  On  the  image,  Muifa  was  a  tropieal  storm  near  the 
eentral  Philippines  at  1611UTC  14  November.  At  this  time,  the  maximum  sustained 
wind  was  estimated  as  18  m/s  (35  knots)  by  JTWC.  The  eloud-band  width  of  TY  Muifa 
was  about  1000  km. 

2.  Entering  the  SCS 

TY  Muifa  was  weakened  further  to  48.9  m/s  (95  knots)  at  OOOOUTC  19 
November  and  slowly  headed  towards  southwest.  It  moved  south-southwestward  slowly 
aeross  Luzon  and  entered  the  SCS.  Its  intensity  was  weakened  to  30.9  m/s  (60  knots)  and 
downgraded  to  a  tropieal  storm  at  0600UTC  20  November. 

At  OOOOUTC  21  November,  TS  Muifa  was  intensified  to  33.4  m/s  (65  knots)  and 
upgraded  baek  to  typhoon.  Further  strengthening  oeeurred  as  TY  Muifa  went  west- 
southwestwards  aeross  the  warm  waters  of  the  SCS.  The  intensity  of  TY  Muifa  reached 
the  seeond  peak  46.3  m/s  (90  knots)  at  1800UTC  at  (11. 6° A,  114.4°E’),  approximately 
800  km  east  of  Vietnam.  TY  Muifa  eontinued  its  journey  to  Vietnam  from  22  to  23 
November  and  was  weakened  again.  The  maximum  veloeity  was  further  deereased  to 
23.2  m/s  (45  knots)  at  1200UTC  24  November. 

3.  Weakening 

Muifa  was  weakened  to  15.4  m/s  (30  knots)  and  downgraded  to  tropieal 
depression  status  at  1200UTC  25  November.  At  OOOOUTC  26  November  Muifa  turned 
northward  into  an  environment  of  inereased  wind  shear  and  as  the  intensity  had  fallen  to 
12.9  m/s  (25  knots),  JTWC  issued  the  final  warning  on  TC  Muifa.  The  final  position  was 
250  km  south-southwest  of  Bangkok,  Thailand. 
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Figure  4  Monthly  tropical  cyclone  frequency  in  the  Western  North  Pacific  Ocean.  The 
upper  line  represents  the  maximum  wind  speed  17  m/sec,  while  the  lower  shadow 
is  the  maximum  wind  Speed  33  m/sec  (from  Neumann  1993). 
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Figure  5  The  best  traek  passage  of  TY  Muifa  (2004)  (from  JTWC  2005). 


Figure  6  The  TRMM  microwave  image  of  TS  Muifa  (2004)  at  161  lUTC  14  November 

close  to  the  central  Philippines  (from  JTWC  2005). 
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Date 

Time 

Lat  (N) 

Long  (E) 

Type 

Pressure 
{/uPa ) 

V 

max 

(m/sec) 

R18  (km) 

R26  (km) 

R33  (km) 

11/16/2004 

00:00 

14.5 

125.7 

TS 

994 

20.6 

148 

— 

— 

11/16/2004 

06:00 

14.5 

124.9 

TS 

991 

23.1 

183 

— 

— 

11/16/2004 

12:00 

14.5 

124.2 

TS 

984 

28.3 

183 

56 

— 

11/16/2004 

18:00 

14.4 

123.6 

TS 

984 

28.3 

174 

56 

— 

11/17/2004 

00:00 

14.6 

123.6 

TS 

980 

30.9 

219 

70 

— 

11/17/2004 

06:00 

14.8 

123.6 

TY 

976 

33.4 

213 

74 

28 

11/17/2004 

12:00 

15.2 

123.8 

TY 

967 

38.6 

213 

74 

28 

11/17/2004 

18:00 

15.5 

123.8 

TY 

954 

46.3 

174 

74 

46 

11/18/2004 

00:00 

15.7 

123.8 

TY 

938 

54 

156 

65 

46 

11/18/2004 

06:00 

15.9 

123.9 

TY 

927 

59.2 

167 

65 

46 

11/18/2004 

12:00 

15.9 

124.2 

TY 

927 

59.2 

170 

70 

46 

11/18/2004 

18:00 

15.7 

124.4 

TY 

938 

54 

133 

59 

41 

11/19/2004 

00:00 

15.2 

124.2 

TY 

938 

54 

109 

59 

41 

11/19/2004 

06:00 

14.7 

124.1 

TY 

954 

46.3 

139 

59 

41 

11/19/2004 

12:00 

14.2 

123.7 

TY 

958 

43.7 

144 

59 

41 

11/19/2004 

18:00 

13.7 

122.8 

TY 

963 

41.2 

131 

56 

37 

1 1/20/2004 

00:00 

12.8 

121.6 

TY 

967 

38.6 

137 

56 

37 

1 1/20/2004 

06:00 

12.5 

120.3 

TY 

976 

33.4 

137 

56 

— 

1 1/20/2004 

12:00 

12.3 

119.3 

TS 

980 

30.9 

115 

56 

— 

1 1/20/2004 

18:00 

12.2 

118.3 

TY 

976 

33.4 

115 

56 

— 

11/21/2004 

00:00 

11.9 

117.2 

TY 

976 

33.4 

128 

56 

— 

11/21/2004 

06:00 

11.9 

116.1 

TY 

972 

36 

152 

56 

28 

11/21/2004 

12:00 

11.8 

115.2 

TY 

963 

41.2 

152 

65 

37 

11/21/2004 

18:00 

11.6 

114.4 

TY 

954 

46.3 

152 

65 

37 

1 1/22/2004 

00:00 

11.4 

113.6 

TY 

958 

43.7 

152 

65 

37 

1 1/22/2004 

06:00 

11.1 

113.1 

TY 

958 

43.7 

126 

59 

33 

1 1/22/2004 

12:00 

10.8 

112.6 

TY 

967 

38.6 

120 

52 

33 

1 1/22/2004 

18:00 

10.5 

112.1 

TY 

976 

33.4 

120 

52 

33 

1 1/23/2004 

00:00 

10.1 

111.7 

TY 

976 

33.4 

120 

52 

— 

1 1/23/2004 

06:00 

9.9 

111.1 

TY 

976 

33.4 

111 

56 

— 

1 1/23/2004 

12:00 

9.6 

110.6 

TS 

980 

30.9 

109 

48 

— 

1 1/23/2004 

18:00 

9.3 

110.2 

TS 

984 

28.3 

107 

44 

— 

1 1/24/2004 

00:00 

9.1 

109.7 

TS 

984 

28.3 

107 

44 

- 

1 1/24/2004 

06:00 

8.8 

108.8 

TS 

987 

25.7 

120 

52 

— 

1 1/24/2004 

12:00 

8.5 

107.4 

TS 

994 

20.6 

102 

- 

- 

1 1/24/2004 

18:00 

8.3 

105.7 

TS 

994 

20.6 

93 

- 

- 

1 1/25/2004 

00:00 

8.7 

103.6 

TS 

994 

20.6 

93 

— 

— 

1 1/25/2004 

06:00 

8.7 

101.7 

TS 

997 

18 

96 

- 

- 

Table  1  The  best  track  record  of  TY  Muifa  (2004)  from  OOOOUTC  16  November  to 

0600UTC  25  November  (from  JTWC  2005). 
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IV.  TYPHOON  WINDS 


During  the  typhoon  seasons,  the  in-situ  measurements  are  diffieult  to  eonduet. 
The  remotely  sensed  data  [e.g.,  QuikSCAT  (QSCAT)]  may  provide  surfaee  winds. 
However,  due  to  the  spatial  eoverage,  the  strong  rotational  motion  from  TY  Muifa  may 
not  be  represented  in  the  QSCAT  data  (Figure  7).  To  overeome  this  defieieney,  a 
Tropieal  Cyelone  Wind  Profile  Model  (TCWPM)  (Carr  and  Elsberry  1997)  is  used  to 
produee  the  high  resolution  gridded  surfaee  wind  field  for  TY  Muifa. 

A,  TROPICAL  CYCLONE  WIND  PROFILE  MODEL 

TCWPM  eomputes  the  wind  field  for  a  tropieal  eyelone  (TC)  using  a  weighted 
formula  by  averaging  inside  and  outside  of  the  tropieal  eyelone  (Figure  8) 

V  =  (1-^)(V+V)+^V^,  (4.1) 


where  V  is  the  total  wind  field,  =  (m,v)  is  the  wind  veetor  produeed  by  the  TC,  V  is 
the  wind  vector  due  to  the  TC  translation,  and  is  the  background  wind  field.  Here 
(m,v)  are  the  TC’s  radial  and  tangential  velocity  components.  The  weight  s  is  defined 
as  (Chu  et  al.  2000) 


s= 


0.97?o  ’ 


(4.2) 


where  r  is  the  distance  to  the  center  of  the  TC  center,  and  7?^  is  the  effective  radius 
defined  by 

v(7?o)  =  0.  (4.3) 

1.  Tangential  Wind  Distribution 

Using  the  angular  momentum  conservation,  the  absolute  angular  momentum 
( )  of  a  uniform  circular  vortex  is  given  as 


Ma{r,p,t)  =  rv  +  -f^^r  , 


2' 


(4.4) 
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where  /g  is  the  Coriolis  parameter  with  respeet  to  the  latitude  of  the  TC  eenter.  In  a 

symmetrieal  TC,  the  absolute  angular  momentum  of  an  air  pareel  depends  only  on  the 
frietional  torque 


dt 


=  rF, 


0’ 


(4.5) 


where  Fg  is  the  frietional  effeet  in  the  tangential  direetion  {0).  In  the  low-level  inflow 
region  of  the  TC,  the  effeet  of  the  frietional  torque  deereases  the  angular  momentum  over 
time  while  it  spiraling  inward  to  the  center  of  the  TC. 


For  slow  temporal  change  of  the  TC  circulation,  depends  on  r  and  pressure 
ip)  only.  Eq.  (4.4)  may  be  rewritten  as 


v{r,p) 


M,{r,p) 

r 


(4.6) 


Since  the  low-level  decreases  faster  near  the  center  and  slower  away  from  the  center, 
is  approximately  represented  by 

MSr,p)  =  M{py-\  (4.7) 


where  X  is  a  positive  constant,  Carr  and  Elsberry  (1997)  proposed  to  use  X  =  0.4  . 
Substituting  (4.7)  into  (4.6)  gives 


v(r,p) 


M{p)  1  r 

2 ■ 


(4.8) 


To  compute  the  tangential  velocity  v{r,p),  the  value  of  M{p)  must  be  determined. 
Substituting  r  =  y  into  (4.8)  and  using  (4.3),  we  have 


(4.9) 

Substituting  (4.9)  into  (4.8)  and  applying  the  adjustment  factora ,  the  tangential  wind  is 
given  by 
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vir.R,.R,)  =  ^ 


(4.10) 


a 


4 


1-a 


4  ’ 


where  is  the  seale  radius  from  the  center  of  the  TC. 

The  radial  wind  of  the  TC  is  computed  from  the  tangential  wind 

M(r)  =  cot(f  )v(r),  (4.11) 

where  /  is  the  inflow  angle  of  the  wind  as  it  spirals  into  the  center  of  the  cyclone  (Figure 
8).  Schwerdt  et  al.  (1979)  suggested  a  relationship  between  the  wind  inflow  angle  /  and 
the  ratio  of  the  radius  to  the  maximum  wind  radius  (7?„,)  (Figure  9). 

2.  Determination  of  Model  Parameters 

In  the  TCWPM,  as  described  in  (4.10),  the  tangential  wind  depends  on  the  radial 
distance  (r),  the  scale  radius  (7?,),  and  the  effective  radius  (7?o)-  Carr  and  Elsberry 
(1997)  determine  from  the  satellite  image  according  to  the  size  of  overall  convective 

and  outflow  cloud  pattern  of  the  TC.  On  the  other  hand,  Chu  et  al.  (2000)  fixed  the 
values  of  radii  for  the  whole  typhoon  passage.  Since  there  is  no  available  satellite  image 
of  TY  Muifa  in  the  SCS,  the  values  of  radii  of  TY  Muifa  are  determined  using  the  data 
from  JTWC  (Table  1).  Near  the  typhoon  center,  the  winds  reported  by  JTWC  can  be 
approximately  taken  as  the  tangential  winds.  Let  tangential  winds  at  (rj,r2,...,r^)  be 

(vj,V2,...,v^.)  .  Eq.  (4.10)  gives  a  set  of  nonlinear  algebraic  equations 


v{r2,Ro,RJ  =  V2, 


(4.12) 


An  optimization  scheme  is  used  to  determine  the  values  of  radii  (R^,R^),  and  to 
further  compute  .  Eor  example,  {R^,Rq)  are  estimated  as  (23.7  km,  795  km)  at 
0600UTC  24  November.  Table  2  shows  temporally  varying  (R^^,R^,Rq)  .  In  this  table, 

the  translation  speed  of  TY  Muifa  is  also  calculated  from  the  distance  between  typhoon 
center  positions  reported  by  JTWC. 
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B,  TCWPM  APPLICATION  ON  TY  MUIFA  (2004) 

1.  QuikSCAT  Satellite 

The  QSCAT  satellite  was  launched  on  19  June  1999  by  NASA  with  the  SeaWinds 
scatterometer  onboard.  QSCAT  has  a  sun-synchronous  orbit  at  an  altitude  of  803  km  and 
a  period  of  101  minutes.  The  ascending  equator  crossing  time  of  it  is  0600UTC. 
SeaWinds  consists  of  a  1  m  long,  18  revolutions  per  minute  rotating  parabolic  antenna, 
with  two  offset  feeds  that  generate  two  13.4-GHz  pencil  beams  at  different  incidence 
angle.  Its  rotating  dish  antenna  sends  out  a  pair  of  scanned  beams  and  ranges  bins  the 
return  from  each  beam  (Martin  2004).  The  detail  description  about  the  design  and 
operation  of  QSCAT  and  SeaWinds  could  be  referred  to  Spencer  et  al.  (1997,  2000)  and 
Liu  (2002). 

The  QSCAT  Level  3  dataset,  which  consists  of  grid  values  of  scalar  and 
componential  wind  velocity,  is  provided  by  the  Jet  Propulsion  Laboratory  (JPL).  The 
Level  3  data  is  produced  on  an  approximately  1/4°  x  1/4°  global  grid  with  both  ascending 
and  descending  passes.  Data  is  available  online  in  Hierarchical  Data  Format  (HDF)  form 
19  July  1999  to  present.  The  QSCAT  Level  3  data  were  downloaded  during  the  period  of 
TY  Muifa  (2004).  The  evolution  of  QSCAT  wind  field  from  17  to  25  November  2004 
(Figure  7)  clearly  shows  that  the  typhoon  strength  agrees  properly  with  the  progressing  of 
TY  Muifa  into  SCS,  but  the  structure  of  tropical  cyclone  winds  is  not  presented  well  due 
to  temporal  limitation  of  satellite  tracks  and  inaccuracy  of  scatterometer  in  high  winds. 

2.  Combination  of  QSCAT  and  TCWPM  Winds 

When  the  background  wind  vector  term  (V^^)  in  (4.10)  is  taken  as  the  temporally 

averaged  QSCAT  Level  3  winds  from  both  ascending  and  descending  passes  during  16- 
25  November  2004  (Figure  10).  The  winds  blow  from  northeast  to  southwest  with  the 
spatial  average  of  about  7.8  m/s.  It  represents  the  dominant  winter  monsoon  during  the 
period  of  TY  Muifa  passage  in  the  SCS. 

With  the  mean  QSCAT  wind  data  as  ,  the  total  wind  field  is  computed  from 

OOOOUTC  16  November  to  0600UTC  25  November  for  (0°  -25  °  A ,  105° A -122° A)  using 
(4. 1 0)  on  a  1  /  4°  x  1  /  4°  grid  with  the  time  interval  of  6  hr.  Such  a  wind  field  is  referred  to 
the  QTCWPM  winds.  The  daily  evolution  of  the  QTCWPM  winds  (Figure  11)  shows 
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that  the  wind  speed  inereases  as  TY  Muifa  enters  the  SCS,  and  deereases  as  Muifa 
approaehes  the  land.  The  QTCWPM  maximum  wind  speed  is  eomparable  to  the 
maximum  wind  speed  reported  in  the  best  traek  reeord  (46.3  m/s).  Furthermore,  the 
QTCWPM  wind  field  shows  asymmetrie  with  the  higher  wind  speed  in  the  right  side. 

The  QTCWPM  winds  are  eompared  to  the  QSCAT  winds  on  1/4°  x  1/4°  as  well 
as  to  the  National  Centers  for  Environmental  Predietion  (NCEP)  surfaee  winds  on  a 
TxT  grid  (Eigure  12).  The  root  mean  square  error  (RMSE)  is  3.2  m/s  between 
QTCWPM  and  NCEP  winds  and  3.4  m/s  between  QTCWPM  and  QSCAT  winds  (Table 
3).  During  TY  Muifa  passage  in  the  SCS,  the  averaged  wind  speed  is  around  35  m/s,  the 
eomputed  QTCWPM  wind  field  is  reasonably  well. 

3,  Ideal  Typhoon  Winds  without  Monsoon 

As  mentioned  above,  the  SCS  baekground  wind  field  during  TY  Muifa  passage  is 
the  winter  monsoon.  In  order  to  separate  and  eompare  typhoon  and  monsoon  foreing  for 
the  oeean  waves,  a  set  of  ideal  typhoon  winds  is  generated  using  (4.10)  with  zero 
baekground  winds  (V^^  =0).  Eigure  13  shows  the  daily  evolution  of  the  asymmetrie 

ideal  typhoon  winds.  The  maximum  wind  speed  is  about  46.3  rn/see.  Before  arrival  and 
after  departure  of  TY  Muifa,  the  wind  speed  is  elosed  to  zero.  Besides,  the  wind  speed  is 
also  zero  far  from  the  typhoon  traek. 
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(a)  17.NOV-2004 


(b)  18-NOV-2004 


(c)  19-NOV-2004 
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Figure  7  Daily  wind  field  of  QSCAT  satellite  observation  on  17-25  November  2004. 
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Figure  8  Wind  field  deeomposition  of  the  TCWPM  model. 


Figure  9  Relation  between  the  wind  inflow  angle  and  the  relative  radius  of  tropical  cyclone 

(after  Schwerdt  1979). 
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Figure  10  Mean  wind  field  of  QSCAT  satellite  observation  from  OOOOUTC  16  November  to 

0600UTC  25  November  2004. 
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(a)  17-NOV-2004 


(b)  18-NOV-2004 


(C)  19-NOV-2004 


(d)  20-NOV-2004 


(e)  21 -Nov-2004 


(f)  22-NOV-2004 


(g)  23-NOV-2004 


(h)  24-NOV-2004 


(i)  25-NOV-2004 


Figure  1 1  Daily  wind  field  of  QTCWPM  calculation  on  17-25  November  2004. 
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Figure  12  Daily  wind  field  of  NCEP  Reanalysis  dataset  on  17-25  November  2004. 
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(b)  18-NOV-2004 


(C)  19-NOV-2004 


(d)  20-NOV-2004 


(g)  23-NOV-2004 


(e)  21 -Nov-2004 


(h)  24-NOV-2004 


(i)  25-NOV-2004 


Figure  13  Daily  wind  field  of  ideal  typhoon  on  17-25  November  2004. 
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Date 

Time 

Position 

TS 

(m/s) 

V 

max 

(m/s) 

Wind  Profile  Radii 

Estimated  Radii 

Lat  (N) 

Long  (E) 

Ris  (km) 

R26  (km) 

R33  (km) 

Rs  (km) 

^max  (km) 

Ro  (km) 

11/16 

00:00 

14.5 

125.7 

3.99 

20.6 

148 

— 

— 

54.8 

89.3 

630 

11/16 

06:00 

14.5 

124.9 

3.99 

23.1 

183 

— 

— 

54.2 

88.9 

680 

11/16 

12:00 

14.5 

124.2 

3.49 

28.3 

183 

56 

— 

43.2 

71.8 

678 

11/16 

18:00 

14.4 

123.6 

3.04 

28.3 

174 

56 

— 

32.6 

54.8 

666 

11/17 

00:00 

14.6 

123.6 

1.03 

30.9 

219 

70 

— 

36.6 

61.5 

725 

11/17 

06:00 

14.8 

123.6 

1.03 

33.4 

213 

74 

28 

18.9 

32.2 

623 

11/17 

12:00 

15.2 

123.8 

2.29 

38.6 

213 

74 

28 

23.8 

40.5 

725 

11/17 

18:00 

15.5 

123.8 

1.54 

46.3 

174 

74 

46 

10.4 

17.9 

638 

11/18 

00:00 

15.7 

123.8 

1.03 

54.0 

156 

65 

46 

8.1 

14.0 

655 

11/18 

06:00 

15.9 

123.9 

1.14 

59.2 

167 

65 

46 

8.1 

14.0 

693 

11/18 

12:00 

15.9 

124.2 

1.49 

59.2 

170 

70 

46 

8.1 

14.0 

693 

11/18 

18:00 

15.7 

124.4 

1.43 

54.0 

133 

59 

41 

8.1 

14.0 

655 

11/19 

00:00 

15.2 

124.2 

2.76 

54.0 

109 

59 

41 

8.1 

14.0 

670 

11/19 

06:00 

14.7 

124.1 

2.62 

46.3 

139 

59 

41 

8.1 

14.0 

615 

11/19 

12:00 

14.2 

123.7 

3.26 

43.7 

144 

59 

41 

9.3 

16.0 

631 

11/19 

18:00 

13.7 

122.8 

5.18 

41.2 

131 

56 

37 

9.8 

16.9 

629 

11/20 

00:00 

12.8 

121.6 

7.59 

38.6 

137 

56 

37 

10.6 

18.2 

661 

11/20 

06:00 

12.5 

120.3 

6.71 

33.4 

137 

56 

— 

15.7 

26.9 

663 

11/20 

12:00 

12.3 

119.3 

5.13 

30.9 

115 

56 

— 

16.8 

28.7 

646 

11/20 

18:00 

12.2 

118.3 

5.06 

33.4 

115 

56 

— 

13.7 

23.5 

648 

11/21 

00:00 

11.9 

117.2 

5.75 

33.4 

128 

56 

— 

14.2 

24.4 

667 

11/21 

06:00 

11.9 

116.1 

5.54 

36.0 

152 

56 

28 

14.7 

25.2 

710 

11/21 

12:00 

11.8 

115.2 

4.56 

41.2 

152 

65 

37 

10.6 

18.3 

714 

11/21 

18:00 

11.6 

114.4 

4.16 

46.3 

152 

65 

37 

8.1 

14.0 

725 

11/22 

00:00 

11.4 

113.6 

4.16 

43.7 

152 

65 

37 

9.1 

15.7 

729 

11/22 

06:00 

11.1 

113.1 

2.96 

43.7 

126 

59 

33 

8.1 

14.0 

718 

11/22 

12:00 

10.8 

112.6 

2.96 

38.6 

120 

52 

33 

10.4 

17.9 

721 

11/22 

18:00 

10.5 

112.1 

2.96 

33.4 

120 

52 

33 

16.7 

28.6 

763 

11/23 

00:00 

10.1 

111.7 

2.89 

33.4 

120 

52 

— 

12.8 

22.0 

726 

11/23 

06:00 

9.9 

111.1 

3.21 

33.4 

111 

56 

— 

11.7 

20.1 

716 

11/23 

12:00 

9.6 

110.6 

2.97 

30.9 

109 

48 

— 

14.3 

24.6 

733 

11/23 

18:00 

9.3 

110.2 

2.55 

28.3 

107 

44 

— 

16.7 

28.6 

738 

11/24 

00:00 

9.1 

109.7 

2.74 

28.3 

107 

44 

— 

16.7 

28.6 

749 

11/24 

06:00 

8.8 

108.8 

4.83 

25.7 

120 

52 

— 

23.7 

40.4 

795 

11/24 

12:00 

8.5 

107.4 

7.29 

20.6 

102 

- 

- 

33.7 

56.9 

775 

11/24 

18:00 

8.3 

105.7 

8.72 

20.6 

93 

— 

— 

30.4 

51.5 

763 

11/25 

00:00 

8.7 

103.6 

10.89 

20.6 

93 

- 

- 

30.5 

51.6 

739 

11/25 

06:00 

8.7 

101.7 

9.67 

18.0 

96 

- 

- 

47.3 

78.8 

772 

Table  2  The  translation  speed  (derived  from  the  best  traek  reeord  as  Table  1)  and  the 

estimated  radii  of  TY  Muifa  (2004). 
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Compared  Datasets 

Wind  Speed 

Wind  Direetion 

Bias  (m/see) 

RMSE  (m/see) 

Bias  (degree) 

RMSE  (degree) 

QTCWPM  /  QSCAT 

-1.902 

3.417 

-2.446 

41.858 

QTCWPM  /  NCEP 

-0.190 

3.228 

-0.287 

34.965 

Table  3  The  statisties  of  QTCWPM  wind  field  eompared  with  QSCAT  observation  and 

NCEP  Reanalysis  dataset. 
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V.  WAVEWATCH-III  MODEL 


A.  MODEL  DESCRIPTION 

WW3  is  a  fully  spectral  third-generation  ocean  wind-wave  model.  It  has  been 
developed  at  the  Ocean  Modeling  Branch  of  the  Environmental  Modeling  Center  of 
NCEP  for  global  and  regional  sea  wave  prediction.  WW3  was  built  on  the  base  of 
WAVEWATCH-I  and  WAVEWATCH-II,  which  were  developed  at  the  Delft  University 
of  Technology  and  NASA  Goddard  Space  Elight  Center,  respectively  (Tolmanl999). 
The  WW3  has  been  validated  over  global-scale  (Wittmann  2001;  Tolman  et  al.  2002)  and 
regional  wave  forecasts  (Chu  et  ah,  2004). 

B,  MODEL  EQUATIONS 

The  governing  equations  and  numerical  approaches  used  in  WW3  are  described 
here.  Eurther  detail  can  be  referred  to  the  WW3  user  manual  (Tolman  1999).  All  the 
notations  used  here  are  following  the  WW3  user  manual.  The  basic  assumptions  in 
WW3  are  slowly  varying  depths  compared  to  an  individual  wave  and  absent  of  currents. 
The  former  assumption  implies  a  large-scale  bathymetry,  for  which  wave  diffraction  can 
generally  be  ignored.  The  latter  suggests  the  energy  of  a  wave  packet  is  conserved 
(Tolman  1999). 

1.  Governing  Equation 

Assuming  incompressible  and  irrotational  flow  at  the  sea  surface,  the  potential 
function  cl)(x,  z,  t)  is  in  a  simplified  relationship  that 

V"O(x,z,0  =  0,  (5.1) 

which  is  also  called  the  Eaplace’s  equation.  Eurthermore,  it  has  a  plane  wave  solution 
that 

0(x,  z.t)  =  sin( Ax  -  crt)/(z),  (5.2) 

/(z)  = ;/ cosh  [A:(z  +  J)] ,  (5.3) 

where  k  is  the  wavenumber,  cr  is  the  relative  radian  frequency,  d  is  the  mean  water 
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depth,  and  ;;  is  a  scale  constant.  The  quasi-uniform  wave  theory  can  be  applied  locally. 
In  a  reference  frame,  the  dispersion  relation  of  the  relative  radian  frequency  {cr  =  2nf^) 
is  given  as 

cr^  =  gk  tanh  kd,  (5 .4) 

(U  =  cr  +  k  •  U,  (5.5) 


where  is  the  reference  frequency,  co  is  the  absolute  radian  frequency,  k  is  the 
wavenumber  vector,  and  U  is  the  averaged  (in  depth  and  time)  current  velocity  vector. 

The  WW3  uses  a  wavenumber-direction  spectrum  F{k,6)  as  a  basic  quantity, 
because  of  its  invariance  characteristics  with  respect  to  physics  of  wave  growth  and 
decay  for  variable  water  depths.  In  the  output,  however,  the  WW3  consists  of  the 
traditional  frequency-direction  spectrum  F(/., 6*)  .  The  transformation  between  two 
spectra  is  calculated  by  the  Jacobean  transformations  (Tolman  1999). 

Because  of  the  wave  action  A  =  E !  <j  is  conserved,  the  WW3  uses  the  wave 
action  density  spectrum  N{k,6)  =  F{k,6)la  within  the  model.  The  wave  propagation 
then  is  described  as 


DN  _S 
Dt  <j  ’ 


(5.6) 


where  D !  Dt  represents  the  total  derivative  of  wave  action  density  spectrum,  and  S 
represents  the  net  effect  of  source  and  sinks  for  the  spectrum  N{k,6) .  This  is  called  the 
balance  equation. 

2,  Wave  Propagation 

In  a  large-scale  application  as  the  SCS,  the  balance  equation  is  transferred  to  a 
spherical  grid,  defined  by  longitude  X  and  latitude  (j)  as 


dN 


dt  cos^ 


\  d  ■  d  ■  d  ■  d  ■  9 

^Ncose  +  —  XN  +  —  kN  +  —  eN  =  -, 


dx 


dk 


dO 


a 


(5.7) 


c^cosO  +  U^ 
R 


(5.8) 
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(5.9) 


sin  6 +  U, 
A  =  ^ ^ 
i?cos^zi 


Cg  tan^zicos6' 


e=o- 

^  R 


(5.10) 


where  R  is  the  radius  of  the  earth,  and  are  eurrent  eomponents.  Eq.  (5.7) 

ineludes  a  eorreetion  term  for  propagation  along  great  eireles. 

3.  Source  Terms 

The  net  source  term  S  is  considered  to  consist  of  four  portions,  a  wind-wave 
interaction  term  5*.^ ,  a  nonlinear  wave-wave  interaction  term  5'^, ,  a  dissipation  (or  called 

white-capping)  term  ,  and  a  additional  wave-bottom  interaction  term  for  shallow 
water.  These  terms  used  in  WW3  are  as 


S  =  S,. 


(5.11) 


These  source  terms  are  defined  for  the  energy  spectra.  Since  in  the  WW3,  most 
source  terms  are  directly  calculated  for  the  action  spectrum,  they  are  denoted  as  S  =  S  /  <j 
in  latter  content.  The  source  term  packages  used  in  this  study  are,  the  Tolman  and 
Chalikov  (1996)  for  both  input  and  dissipation  terms,  the  discrete  interaction 
approximation  (DIA)  method  for  nonlinear  interaction,  and  the  Joint  North  Sea  Wave 
Project  (JONS WAP)  formulation  for  bottom  friction. 

a.  Tolman  and  Chalikov  (1996)  Input  Term 

The  source  terms  package  of  Tolman  and  Chalikov  (1996)  include  input 
and  dissipation  terms.  The  input  source  term  is  given  as 


Si„{k,0)  =  <JpN{k,6), 


(5.12) 


where  is  a  nondimensional  wind-wave  interaction  parameter,  and  can  be  computed  as 


^2’ 

-1  <  <7^  <  Qj  /  2, 

Qj  /  2  <  <  Qj 

a-jG a  ~ 

Qj  <  <7^  <  Qj’ 

(5.13) 


with 
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(J, 


(5.14) 


CTti  1  ,  y-.  y-i  \ 

= - cos(6»-6»„), 


where  is  the  nondimensional  frequeney  of  a  speetral  eomponent,  0^  is  the  wind 
direetion,  and  is  the  wind  veloeity  at  a  height  equal  to  the  apparent  wave  length 


=■ 


Itt 


k  cosi0-0J 


(5.15) 


The  parameters  to  ,  Qj  and  Qj  (5.13)  depend  on  the  drag 
eoeffieient  Q  at  the  height  of  apparent  wave  length 


Qi  =1.075 +  75Q 

Qj  =1.2  +  300Q, 

=  0.25a5  / a^, 

=  0.25  +  395Q, 

02  =  0.35  +  150Q, 

Oj  =  (Oq  -02  -o^)/(Og  -O4  +0^), 

=  0.30  +  300(7^ , 

0^  =o^Q^, 

ag  =ao(l-a3), 

O^  (i3g  (^2  1^  ^10^  ^(^2  ^1 )’ 

Og  =0.35  +  240Q, 

ajo  =-0.05  +  470Q. 


The  wave  model  take  the  wind  at  a  given  referenee  height  (in  this 

study  uses  10  m)  as  its  input,  so  that  and  Q  need  to  be  derived  as  part  of  the 

parameterization.  Exeluding  a  thin  surfaee  layer  adjusting  to  the  water  surfaee  the  mean 
wind  profile  is  elose  to  logarithmie 


u,=^ln(%  (5.17) 

where  ic  =  0.4  is  the  Von  Karman  eonstant,  and  z^  is  the  roughness  parameter.  This 
equation  ean  be  rewritten  in  terms  of  the  drag  eoeffieient  at  the  referenee  height  z^  as 


with 


C^=/c^[R-HC)f, 


R  =  ln 


^rg 


X-Jaul 


(5.18) 


(5.19) 
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where  x  -  0-2  is  a  eonstant,  and  a  is  the  eonventional  nondimensional  energy  level  at 
high  frequeneies.  An  aecurate  explieit  approximation  to  these  implieit  relations  is  given 
as 


c,  =  io-' 


( 

0.021  + 

V 


10.4  ^ 


(5.20) 


The  estimation  of  the  drag  eoeffieient  thus  requires  an  estimate  of  the 
high-frequeney  energy  level  a  ,  whieh  could  be  estimated  directly  from  the  wave  model. 
However,  the  corresponding  part  of  the  spectrum  is  generally  not  well  resolved,  tends  to 
noisy,  and  it  tainted  by  errors  in  several  source  terms.  Therefore,  a  is  estimated 
parametrically  as  (Janssen  1989) 


a  =  0.57 


7 

tl* 


\^P  J 


(5.21) 


Using  the  definition  of  drag  coefficient  and  (5.20),  the  roughness  parameter  Zq  becomes 

z  o=z^exp(-x-C;‘^^),  (5.22) 


and  the  wind  velocity  and  drag  coefficient  at  height  become 

ln(;L^/zo) 

M,  =U^ - - - 

ln(A/z„) 


C  =C 

'-'A 


r 


\^A  ) 


(5.23) 


(5.24) 


Finally,  (5.21)  requires  an  estimate  for  the  peak  frequency  .  To  obtain  a  consistent 

estimate  of  the  peak  frequency  of  actively  generated  waves,  even  in  complex  multi-modal 
spectra,  this  frequency  is  estimated  from  the  equivalent  peak  frequency  of  the  positive 
part  of  the  input  source  term 

msLx\0,Si^{k,6)\dfd6 
jj  max  [0,  {k,  6)]  dfdO  ’ 

from  which  the  actual  peak  frequency  is  estimated  as 
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(5.25) 


4  =  3.6 X  lO  V  0.924,  - 6.3 X 


(5.26) 


Furthermore,  to  reduce  the  swell  due  to  opposing  or  weak  winds,  a  fdtered 
input  source  term  is  defined  as 

S,n  for  or 

S,„=\XA„  for  /3<0  and  0.6f^<f<0.Sf^,  (5.27) 

for  y0<O  and  f<0.6f^, 

where  /  is  the  frequency,  4  is  the  peak  frequency  of  the  wind  sea  as  computed  from 
the  input  source  term,  and  0  <  <  1  is  a  reduction  factor  for  Si„ ,  which  is  applied  to 

swell  with  negative  P .  represents  a  linear  reduction  of  with  4  providing  a 

smooth  transition  between  the  original  and  reduced  input. 

b.  Tolman  and  Chalikov  (1996)  Dissipation  Term 
The  dissipation  source  term  consists  of  two  constituents.  The  dominant 
low-frequency  constituent  is  based  on  an  analogy  with  energy  dissipation  due  to 
turbulence, 

^ds,i  =  -'lu»hk^(/)N{k,  6), 

h  =  4\^l^[j(f,d)dfdd^  ,  (5.28) 

where  h  is  a  mixing  scale  determined  from  the  high-frequency  energy  content  of  the 
wave  field  and  where  (j)  is  an  empirical  function  accounting  for  the  development  stage  of 

the  wave  field.  By  defining  a  minimum  value  for  (j)  at  minimum  frequency  4 ,  mm  ’ 

the  nonlinear  term  has  been  added  to  allow  for  some  control  over  fully  grown  conditions. 
If  is  below  the  linear  curve,  and  are  given  as 


bl  fp,i,mm  ^  ^^min  ^0  )  ’ 


63  =8. 


(5.29) 


If  is  above  the  linear  curve,  b^  and  p  are  given  as 
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(5.30) 


Z  =max{X-0.0025, 

^2  “  f p,i, mm  ^  ^min  “^0  ~^\fb  ’ 

6,= 

^min  -^0  -bJb 

The  above  estimate  of  results  in  d(/)l  dfp^=Q  for  ,  =  /*  •  For  <  /^ ,  the  (/)  is  kept 
constant  that  (/)  =  . 

The  empirical  high-frequency  dissipation  id  defined  as 

(u  Y 

S,.,bik,e)  =  -a,  ^  fa:N{k,e), 

\s) 

B  =  a,  ^  ,  (5.31) 

I  g  ; 

a„=^^tN{k,0)de, 

Cgg 

where  is  Phillips’ (1957)  nondimensional  high-frequency  energy  level  normalized 
with  ,  and  where  through  and  are  empirical  constant. 

The  two  constituents  of  the  dissipation  source  term  are  combined  using  a 
simple  linear  combination,  defined  by  the  frequencies  /  and  /j . 


1  for  f  <fi 

TTT  Zs/<Z 

J\  Jl 

,0  for  f^<f 


To  enhance  the  smoothness  of  the  model  behavior  for  frequencies  near  the  parametric 
cut-off  ,  the  parametric  high-frequency  tail  as 


N{k,  ,e)  =  {\-  cb)n{K  ,  0) + ,  0)  -f- 

V  Ji-l 


(5.33) 
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where  i  is  a  discrete  wavenumber  counter,  m  is  set  to  5,  and  where  ®  is  defined 
similarly  to  ,  ranging  from  0  to  1  between  and  .  The  frequencies  defining  the 

transitions  and  the  length  scale  h  are  predefined  in  the  model  as 


A/ =3.004,, 
/  =1-754,, 
/2  =2.504,, 
fu  =2.004,, 


(5.34) 


c.  Nonlinear  Interaction  Term 

Nonlinear  wave-wave  interactions  are  modeled  using  the  DIA  method 
(Hasselmann  et  al.  1985).  This  parameterization  was  originally  developed  for  the 
spectrum  F{f^,6)  .  To  assure  the  conservation  nature  of  4/  for  this  spectrum,  this 


source  term  is  calculated  for  F{f^,6)  instead  of  N{k,6) . 


Resonant  nonlinear  interactions  occur  between  four  wave  components 
(quadruplets)  with  wavenumber  vector  kj  through  k, .  In  the  DIA,  it  is  assumed  that 

kj  =  kj .  Resonance  conditions  then  require  that 


1^1  +1^2  —  k3  +k4, 
(J2  =  CTj 

CTl  =(i  +  ^M, 

,  0-4  =  (I -4)0-1, 


(5.35) 


where  is  a  constant.  For  these  quadruplets,  the  contribution  dS^i  to  the  interaction 
for  each  discrete  {f^,0)  combination  of  the  spectrum  corresponding  to  kj  is  calculated 
as 


f  2^ 

f 

=  D 

-1 

1 

^  -4  /-ll 

Cs 

4' 

1-4 

V 

4 

(l  +  4z)^ 


+ 


(l-4rj 


444 

(i-4z)^ 


(5.36) 
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where  Fy=  F{f^^,6y)  and  ,  and  so  others,  C  is  a  proportionality 

constant.  In  this  study  the  values  suggested  by  Tolman  and  Chalikov  (1996)  is  applied 
that =  0.25  and  C  =  1.00x1 0\ 

The  water  depth  is  scaled  by  the  factor  D  as 


D  =  \  +  ^ 
kd 


\-c^kd  \e 


kd 


(5.37) 


Recommended  values  for  these  constants  are  Cj=5.5  ,  C2=5/6  ,  and  03=  1.25 

(Hasselmann  and  Hasselmann  1985).  The  over-bar  notation  denotes  straightforward 
averaging  over  the  spectrum.  For  an  arbitrary  parameter  z  the  spectral  average  is  given 
as 

z=E-^r\\F{i,e)dide, 

°  “  (5.38) 

E-Ciy{f,.e)df,d0. 

For  numerical  reasons,  however,  the  mean  relative  depth  is  estimated  as 

kd  =  0.15kd,  (5.39) 


where  k  is  defined  as 


k  =  |l/Vkj".  (5.40) 

d.  JONSWAP  Bottom  Friction  Term 

The  JONSWAP  parameterization  of  bottom  friction  is  simply  represented 
by  (Hasselmann  et  al.  1973) 

S^^^{k,e)  =  2T^^^N{k,e),  (5.41) 

gd 

where  F  is  an  empirical  constant,  which  is  estimated  and  as  F  = -0.061 for  wind 
sea  (Bouws  and  Komen  1983). 
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C.  MODEL  VERIFICATION 

The  WW3  was  integrated  for  the  SCS  from  3  January  to  3 1  Deeember  2000  using 
twiee  daily  gridded  QSCAT  oeean  winds,  and  verified  (Chu  et  al.  2004)  using  the 
signifieant  wave  height  (H^)  obtained  from  T/P  satellite  and  the  European  Remote 

Sensing  Satellites  ERS-1/2  (Eigure  14).  The  WW3  model  errors  stratify  the  Gaussian- 
type  distribution  with  the  RMSE  of  0.48  m,  which  is  comparable  to  the  T/P  altimeter 
accuracy  0.5  m.  This  shows  the  capability  of  WW3  for  SCS  wave  simulation. 
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(a) 


(b) 


5 


Figure  14  Monthly  mean  signifieant  wave  height  (H^)  from  (a)  WW3  simulations,  and  (b) 

T/P  observations  (from  Chu  et  al.  2004). 


39 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


40 


VI.  RELATIONSHIPS  BETWEEN  WAVES  AND  WINDS 


Surface  winds  have  large  variability  in  spaee  and  time.  Different  winds  generate 
different  wave  patterns.  For  slowly  varying  winds,  ehange  of  waves  follows  the  ehange 
of  winds  (loeal  equilibrium).  For  fast  varying  winds  sueh  as  tropical  cyclones,  ehange  of 
waves  does  not  follow  the  ehange  of  winds  (non-loeal  equilibrium).  Wave  field  must 
respond  to  the  direetional  ehange  as  well  as  the  magnitude  variation  of  wind  forcing 
(Jones  and  Toba  2001). 

A.  GENERAL  RESPONSES 

1.  Wind  Speed 

While  steady  winds  transfer  momentum  through  the  sea  surfaee  via  waves,  the 
wave  field  also  plays  an  important  role  in  this  momentum  transfer  through  motion  and 
adjustment.  As  the  adjustment  is  finished,  the  eontinuity  between  winds  and  momentum 
fluxes  would  be  achieved  in  both  air  and  sea  boundary  layers.  Toba  (1972)  referred  this 
balance  eondition  as  the  “loeal  equilibrium.” 

Toba  et  al.  (1988)  reanalyzed  the  dataset  of  Shirahama  (also  in  Kawai  et  al.  1977) 
to  investigate  the  response  of  one-dimensional  wave  frequeney  speetra  to  unsteady  winds. 
Figure  15  shows  the  one-dimensional  wave  speetra  under  increasing  wind  (Figure  15a), 
steady  wind  (Figure  15b),  and  decreasing  wind  (Figure  15c).  The  tendency  of  the  peak 
can  be  clearly  seen  in  all  three  panels.  For  a  rapid  inerease  in  wind  speed  (Figure  15a), 
the  wave  field  becomes  under-saturated.  The  energy  decreases  near  the  spectral  peak. 
Higher  wind  speed  also  inereases  the  amplitude  of  waves,  and  the  higher  frequeney 
waves  inerease  and  remove  momentum  from  the  air.  When  the  reduetion  in  wind  speed 
oeeurs  (Figure  15e),  the  wave  field  is  over-saturated.  The  peak  energy  inereases  and 
beeomes  steeper.  The  high-frequeney  waves  deerease  and  tend  to  extraet  less  momentum 
from  the  air. 

A  further  study  was  conduced  from  the  laboratory  experiments  of  Waseda  et  al. 
(2001).  The  speetral  development  for  a  rapid  ehange  between  low  winds  (4.6  m/s)  and 
high  winds  (7.1  m/s)  is  shown  in  Figure  16.  Two  wave  adjustment  states  exist  for 
increasing  and  deereasing  winds.  During  the  rapid  adjustment  stage,  changes  of  low- 
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frequency  wave  energy  are  not  evident  whereas  the  changes  of  high-frequency  wave 
energy  are  evident.  After  the  rapid  adjustment  period,  the  wave  field  evolves  to  a 
conventional  growth  and  propagation  stage.  As  the  wind  reduces,  the  peak  frequency 
shifts  toward  the  high  frequency  domain;  the  steepness  and  peak  of  the  wave  spectrum 
reduce.  On  the  other  hand,  as  the  wind  enhances,  the  peak  frequency  shifts  toward  the 
low  frequency  domain;  the  steepness  and  peak  of  the  wave  spectrum  enhance. 

2.  Wind  Direction 

The  properties  of  directional  wave  spectra  under  veering  wind  were  investigated 
by  several  measurements  (Hasselmann  et  al.  1980;  Allender  et  al.  1983;  Walsh  et  al. 
1989;  Masson  1990).  The  high-frequency  waves  respond  more  quickly  to  a  directional 
wind  shift  than  the  low-frequency  waves.  Small  wind  shift  gradually  rotates  the  wave 
spectrum  in  accordance  to  the  new  wind  direction.  Large  wind  shift  generates  a  second 
wave  packet  in  the  new  wind  direction  with  gradually  decaying  old  wave  packet.  After 
the  local  equilibrium  between  the  wave  spectra  and  winds  is  established,  the  old  wave 
energy  still  travels  in  the  initial  direction.  This  can  be  treated  as  swell  at  an  angle  to  the 
younger  wind- waves. 

Figure  17  shows  the  directional  evolution  of  five  wave  packets  with  inshore  and 
offshore  wind  shifts  in  an  open  ocean.  The  wave  response  to  the  changing  wind  direction 
strongly  depends  on  the  frequency. 

The  misalignment  of  the  atmospheric  surface  layer  wind,  stress  vector,  and 
dominant  wave  direction  were  also  reported  in  field  observations  (Smith  1980;  Geemaert 
1988;  Long  et  al.  1994).  The  directional  wave  spectra  derived  from  Surface  Contour 
Radar  (SCR)  (Figure  18)  shows  that  the  surface  wind  stress  vector  and  the  dominant 
wave  direction  do  not  always  parallel  the  mean  wind  direction,  and  that  the  wave  spectra 
are  asymmetric  left  and  right  of  the  wind.  Rapid  change  in  wind  direction  causes  initially 
non-local  equilibrium  between  the  wave  spectra  and  winds  (Jones  and  Toba  2001). 

3.  Strong  Cyclonic  Forcing 

Observational  studies  were  conducted  on  the  wave  characteristics  under  strong 
cyclonic  forcing.  Wyatt  (1995)  measured  the  directional  spectra  of  storm  waves  and  to 
examine  the  effect  of  wave  fetch  in  the  Celtic  Sea  using  the  high-frequency  radar.  Holt  et 

al.  (1998)  determined  the  wave  field  from  an  intense  storm  using  the  synthetic  aperture 

42 


radar  (SAR)  imagery  from  ERS-1  satellite.  Wright  et  al.  (2001)  investigated  the  spatial 
variation  of  direetional  wave  speetra  during  hurrieane  Bonnie  using  the  NASA  Scanning 
Radar  Altimeter  (SRA).  They  found  that:  (1)  wind  change  may  cause  rapid  change  in 
short-wave  direction  but  not  in  the  mean  direction  of  the  long-waves;  and  (2)  the  highest 
and  longest  waves  occur  in  the  right  forward  quadrant  of  the  hurricane  (Figure  19). 

B.  NUMERICAL  MODELING 

1.  Source  Terms 

As  mentioned  earlier,  the  source  term  package  of  Tolman  and  Chaltkov  (1996)  is 
chosen  for  the  WW3  modeling  in  this  study.  The  chosen  wind-input  term  {Si„)  has 
several  improved  features, 

•  It  becomes  negative  for  waves  traveling  at  large  angle  with  the  wind  or 
faster  than  the  wind. 

•  It  results  in  a  two  to  three  times  smaller  integral  input  of  energy  for  fully 
grown  sea. 

•  At  high  frequencies  it  results  in  a  larger  energy  input  than  previous 
relation,  and  the  integral  input  for  young  waves  is  more  consisted  with 
previous  relation. 

Also,  the  wave  dissipation  term  ( )  is  divided  into  two  constitutes  according  to 
different  characteristics  in  the  frequency  domain.  In  low-frequency  domain,  is 
calculated  using  the  wave  energy  dissipation  due  to  oceanic  turbulence.  In  high- 
frequency  domain,  is  parameterized  using  the  assumption  of  quasi-steady  balance  of 

source  terms  (Tolman  and  Chaltkov  1996). 

2.  WW3  Model  Simulation 

Recently,  Moon  et  al.  (2003,  2004a,  2004c)  investigated  the  effect  of  hurricane 
wind  forcing  on  surface  waves  and  air-sea  momentum  exchange  using  WW3  with  various 
wind  fields.  These  wind  fields  are  either  idealized  or  computed  from  observational  data 
(during  hurricane  Bonnie  1998). 

Figure  21  shows  the  directional  wave  spectra  at  four  points  on  the  circle  around 
the  storm  center  with  a  radius  of  .  Three  different  hurricane  translation  speeds  (HTS), 

0  m/s,  5  m/s,  and  10  m/s,  were  used  in  the  calculation.  As  the  storm  moves  faster,  the 
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wave  spectra  in  the  east  and  the  north  locations  show  a  unilateral  shape  due  to  resonance 
effect,  while  the  spectra  in  the  west  and  the  south  locations  have  a  more  complex 
structure. 

The  spatial  distribution  of  and  dominant  wave  direction  under  various  HTS  is 

represented  in  Figure  22.  As  HTS  increases,  waves  in  the  front-right  quadrant  of  the 
storm  track  become  higher  and  longer,  while  those  in  the  rear-left  quadrant  become  lower 
and  shorter.  That  is  to  say  the  hurricane-generated  wave  field  is  mostly  determined  by 
two  factors:  the  distance  from  the  hurricane  center  and  the  HTS. 

Figure  23  shows  the  relation  of  the  drag  coefficient  (Q)  with  the  wave  age  and 

the  driving  wind  speed.  It  is  found  that  the  drag  coefficient  depends  on  both  the  wind 
speed  and  the  wave  age;  it  decreases  with  the  wave  age  at  lower  wind  speed  but  increases 
with  the  wave  age  at  very  high  wind  speed  (above  40  m/s)  (Moon  et  al.  2004b). 
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Figure  15  Ensemble  means  of  ten  raw  spectra  under  (a)  increasing  wind  (triangle),  (b) 
steady  wind  (circle),  and  (c)  decreasing  wind  (inverse  triangle)  (from  Toba  et  al. 

1988). 
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Figure  16  Spectral  development  for  sudden  change  of  wind  Speed,  (a)  increasing  wind  case, 
and  (b)  decreasing  wind  case.  Time  steps  are  initial  equilibrium  (solid  line),  1 
second  after  change  (square  mark),  4  second  after  change  (solid  circle  mark),  and 
final  equilibrium  (dashed  line),  respectively  (from  Waseda  et  al.  2001). 


45 


0  S  10  15  20  25 

time  (hour) 

Figure  17  Wind  direction  (solid  line  no  symbol)  and  mean  wave  direction  at  frequency  of 
0.10  (plus  mark),  0.20  (star),  0.31  (diamond),  and  0.54  (square)  Hz,  under  (a) 
inshore  wind,  and  (b)  offshore  wind  (from  Masson  1990). 


Figure  18  Contour  plot  of  directional  wave  spectra.  The  radial  coordinate  is  frequency  in 
Hz,  and  the  arrow  indicates  the  wind  direction  (from  Long  et  al.  1994). 
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Figure  19  Spatial  distribution  of  significant  wave  height  (H^)  measured  by  the  SRA  during 

Hurricane  Bonnie  on  24  August  1998.  At  that  time,  Hurricane  Bonnie  is  moving 
northward  (from  Wright  et  al.  2001). 


Figure  20  Relationship  of  drag  coefficient  ( Q  )  and  sea  surface  wind  speed  (  U^q  )  under 
extreme  conditions.  Different  symbols  indicated  the  Q  derived  from  different 

formulas  (from  Powell  et  al.  2003). 
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Figure  21  Directional  wave  spectra  under  idealized  hurricanes  case  with  different  HTS. 

Locations  are  at  away  from  the  storm  center  in  directions,  (a)  East,  (b)  North, 

(c)  West,  and  (d)  South,  respectively  (from  Moon  et  al.  2004a). 


48 


D»tanca  North  of  Cya  i«*n} 


(a)  Wave  F«lds  lO  m  s  I  (b)  Wavt- FkjWs  15  m  s  ■  (c)  Wave  Felds  lOms 


Ostance  Cast  C)r0  fkT'}  Distance  East  ul  Eye  iKn^  Oisl^r-oe  EsssI  of  Eye  iKm) 


Figure  22  Spatial  distributions  of  significant  wave  height  (H^)  under  various  HTS.  Arrows 


indicate  the  local  mean  wave  direction  (from  Moon  et  al.  2004a). 
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Figure  23  Scatter  plot  of  with  various  wave  age  and  wind  speed.  Q  is  plotted  with 
different  colors  according  to  wind  speed  intervals  (from  Moon  et  al.  2004b). 
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VII.  SCS  WAVE  CHARACTERISTICS  DURING  TV  MUFIA 


The  WW3  model  simulation  of  the  SCS  waves  is  performed  under  three 
seenarios.  First,  the  QTCWPM  wind  field  of  Typhoon  Muifa  (2004),  as  mentioned 
earlier,  is  used  to  drive  the  WW3  model  for  the  SCS.  The  WW3  result  is  evaluated  with 
the  signifieant  wave  height  measured  by  the  TOPEX/Poseidon  satellite.  Seeond,  the  ideal 
typhoon  wind  field  is  also  applied  to  simulate  the  ideal  typhoon-generated  waves  in  the 
SCS.  Third,  the  same  ideal  typhoon  wind  field  with  a  uniform  water  depth  is  run  in 
WW3  to  examine  the  topography  effeet.  The  wave  eharaeteristies  in  the  SCS  is  analyzed 
and  diseussed  by  the  signifieant  wave  height  and  direetional  wave  speetra. 

A.  WW3  MODEL  VALIFICATION  DURING  TY  MUIFA 

1.  Model  Setting  and  Calculation 

This  study  uses  a  relatively  high  resolution  setting  for  the  WW3  model  (Wittmann 
2001;  Tolman  et  al.  2002).  The  spatial  grid  resolution  is  1/4°  x  1/4°  in  the  SCS  of 
latitude  from  0°  to  25°  N ,  and  longitude  from  105  "E  to  122°E’.  The  diseretization  of 
energy  speetra  is  in  24  direetions  with  15  °  grid  spaeing,  and  in  25  frequeneies  with 
logarithmie  inerement  from  0.0418  to  0.3058  Hz.  Two  values  are  used  in  the  WW3  time 
steps  for  eomputational  effieieney.  The  global  time  step,  the  spatial  time  step,  and  the 
speetral  time  step  are  all  set  to  300  seeond,  for  entire  solution  propagation,  spatial 
propagation,  and  intra-speetral  propagation,  respeetively.  The  souree  time  step  for  the 
souree  term  integration  is  set  to  100  seeond. 

The  QUICKEST  with  UETIMATE  propagation  seheme  is  used  here.  As 
previously  diseussed,  four  souree  term  eomponents  are  Tolman  and  Chalikov  (1996) 
wind-wave  interaetion  term  and  dissipation  term,  the  DIA  nonlinear  interaetion  term,  and 
the  JONS  WAP  bottom  frietion  term,  respeetively.  Other  important  model  switeh 
parameters  are  listed  in  Table  4.  The  WW3  model  in  this  study  is  based  on  an  available 
version,  whieh  has  been  implemented  in  the  Naval  Oeean  Analysis  and  Predietion 
Eaboratory  (NOAP)  (Chu  et  al.  2004). 

The  bathymetry  data  is  from  Terrain  Base  (Eigure  24),  whieh  is  the  best  available 
publie  domain  global  model  of  terrain  and  bathymetry  on  a  1/12°  x  1/12°  grid.  This 


51 


dataset,  produced  by  the  National  Geophysical  Data  Center  and  World  Data  Center-A  for 
Solid  Earth  Geophysics,  is  from  twenty  different  sources,  including  the  DCW-DEM  and 
ET0P05.  The  QTCWPM  wind  field  is  used  as  the  input  wind  to  run  WW3  for  the  SCS. 
The  WW3  model  is  run  from  OOOOUTC  16  November  to  1200UTC  25  November  2004. 
Since  TY  Muifa  entered  the  SCS  as  late  as  19  November,  the  model  computation  of  the 
first  three  days,  from  16  to  18  November,  could  be  considered  as  the  ‘spin  up’  period  of 
WW3  model. 

2.  TOPEX/Poseidon  Satellite 

Because  of  the  lack  of  direct-measured  platforms,  such  as  vessels  or  buoys,  during 
TY  Muifa  (2004)  in  the  SCS,  other  sources  of  wind-wave  measurement  must  be 
considered  for  evaluating  WW3.  Satellite  remote  sensing  is  an  important  field 
observation.  Eor  the  remarkable  performance  of  radar  altimetry  and  for  the  data 
availability,  the  altimetry  onboard  the  T/P  satellite  is  applied  in  this  study  (Chu  et  al. 
2004). 

The  T/P  satellite,  designed  to  map  ocean  surface  topography,  is  jointly  launched 
by  NASA  and  the  Trench  Space  Agency,  Centre  National  d’ Etudes  Spatiales  (CNES). 
The  satellite  first  served  in  August  1992,  and  originally  was  planed  for  three-year  prime 
mission,  then  expended  for  another  six-year  extension,  totally  for  more  than  ten  years.  It 
carried  a  dual-frequency  radar  altimeter  and  was  maneuvered  in  a  circular  orbit  at  an 
altitude  of  1336  km  with  a  9.916-day  repeat  period  (Martin  2004). 

The  T/P  altimeter  measurement  data  is  also  available  online.  The  T/P  dataset  is 
provided  by  the  Physical  Oceanography  Distributed  Active  Archive  Center  (PO.DAAC) 
of  JPE.  Data  is  in  merged  GDR  format  generation  B,  product  number  068,  and  organized 
in  on  cycle  header  file  with  maximum  254  pass  files.  The  T/P  dataset  includes  altitude, 
altimeter  range,  and  significant  wave  height  measurements  from  both  Ku-band  and 
C-band. 

3,  Result  Evaluation 

During  the  period  of  TY  Muifa  (2004),  the  temporal  covered  cycles  are  448, 
which  from  11  to  21  November,  and  449,  from  21  November  and  01  December.  There 
are  14  passes  across  the  SCS  area,  which  include  001,  012,  038,  051,  064,  077,  088,  114, 
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127,  140,  153,  164,  216,  and  229.  These  passes  have  total  25  erossover  points  (see 
Figure  25)  in  the  ocean  region.  On  each  point,  the  from  T/P  measurement  is 

accumulated  to  evaluate  the  from  WW3  computation. 

The  comparison  between  WW3  result  and  T/P  observation  is  computed  and 
performed  as  Figure  26.  The  total  number  of  pairs  during  the  period  of  TY  Muifa  is 

38.  The  scattering  diagram  (Figure  26a)  shows  the  clustering  of  points  approximately 
located  around  the  equal  line.  The  histogram  of  the  differences  between  model  and 
observation  (Figure  26b)  shows  a  Gaussian- type  distribution.  Examining  the  statistics, 
the  bias  of  two  datasets  is  0.137  m  with  the  root  mean  square  error  is  0.308  m,  and  the 
correlation  coefficient  is  0.895.  The  statistics  suggests  that  the  WW3  result  and  the  T/P 
observation  agreed  well,  and  the  numerical  simulation  during  the  period  of  TY  Muifa  is 
accurate  and  reasonable. 

B,  NUMERICAL  SIMULATION  SCHEME 

Previous  studies  has  shown  that  the  distance  from  the  typhoon  center,  the 
hurricane  (typhoon)  intensity  and  the  hurricane  (typhoon)  translation  speed  are  the 
determinants  of  the  typhoon-generated  wave  field  (Walsh  et  al  1989;  Wright  et  al.  2001; 
Moon  et  al.  2003).  Several  typhoon  center  locations  are  designed  here  to  investigate  the 
effects  of  typhoon  intensity  and  of  typhoon  translation  speed  on  the  SCS  wave 
characteristics.  Furthermore,  one  set  of  ideal  typhoon  winds,  and  one  set  of  uniform 
water  depth  are  also  used  here  to  investigate  the  effects  of  winter  monsoon  and  of  the 
SCS  topography. 

1.  Typical  Locations  along  the  Typhoon  Track 

Since  the  intensity  and  the  translation  speed  of  TY  Muifa  are  varied  in  time,  four 
cases  are  designed  here  by  fixed  either  variables.  The  typhoon  passage  could  be 
reviewed  in  Chapter  III  (Figure  5).  Figure  27  shows  the  translation  speed  and  the 
intensity  of  TY  Muifa  (data  from  Table  2)  along  the  typhoon  passage.  The  water  depth 
(data  from  Figure  24)  of  the  typhoon  center  is  also  accompanied.  The  location  of  these 
designed  typhoon  centers  is  shown  as  Figure  28. 

Along  the  track  of  the  TY  Muifa  passage,  the  translation  velocity,  low  pressure. 
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intensity,  and  bathymetry  vary  (Figure  27).  Four  loeations  along  the  TY  passage  are 
selected  here  for  the  discussions:  Center-I  (arriving  the  SCS),  Center-II  (strengthening 
and  slowing),  Center-Ill  (weakening),  and  Center-I V  (shelf  breaking),  respectively. 

2.  Designed  Typhoon  Centers 

Center-I  is  chosen  at  (II.9°A^,  II7.2°£’).  Just  after  TY  entering  the  SCS  for  24 
hours,  the  typhoon  center  is  located  at  Center-I  at  OOOOUTC  21  November.  The  typhoon 
translation  speed  slightly  decreased  to  5.8  m/s.  The  typhoon  low  pressure  was  remained 
as  low  as  976  //Pa  with  the  maximum  wind  speed  was  33.4  m/s.  At  this  time,  the 

estimated  maximum  wind  radius  (/?„,)  is  24.4  km.  Center-I  is  for  the  typhoon  arriving 
period. 

Center-II  is  chosen  at  (II.6°A,  II4.4°£’).  TY  Muifa  reached  the  strongest 
intensity  (i.e.,  the  lowest  central  pressure,  954  //Pa  )  at  Center-II  at  I800UTC  21 
November.  The  maximum  wind  speed  was  46.3  m/s,  and  was  14.0  km.  At  that  time, 

the  typhoon  translation  speed  slowed  down  to  4.2  m/s.  Center-II  is  for  the  maximum 
typhoon  intensity  period. 

Center-Ill  is  chosen  at  (10.5  ° N ,  1 12. 1  ° E).  TY  Muifa  had  the  same  intensity  at 
Center-Ill  at  I800UTC  22  November  as  at  Center-I  at  OOOOUTC  21  November.  The 
typhoon  translation  speed  is  notably  slower  at  Center-Ill  (3.0  m/s)  than  at  Center-I  (5.8 
m/s).  The  maximum  velocity  radius  is  slightly  larger  at  Center-Ill  (7?^,  28.6  km)  than  at 

Center-I  24.4  km).  Other  conditions  are  the  same  at  Center-Ill  (minimum  pressure 
976 //Pa,  maximum  wind  speed  33.4  m/s)  as  at  Center-I.  Center-Ill  is  for  the  typhoon 
weakening  period.  Comparison  between  Center-I  and  Center-Ill  leads  to  the  effect  of  the 
typhoon  translation  speed. 

Center-IV  is  chosen  at  (8.8 °A,  108.8 °£’).  TY  Muifa  was  close  to  the  Asian 
landmass  at  0600UTC  24  November.  The  water  depth  is  about  116  m.  The  low  pressure 
was  987  //Pa  ,  and  the  maximum  wind  speed  was  25.7  m/s  with  R^  40.4  km.  The 

typhoon  intensity  is  much  weaker  at  Center-IV  than  at  Center-II.  Center-IV  is  selected 
for  the  typhoon  departure  period. 
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C.  EFFECTS  OF  TYPHOON  WINDS 
1.  Significant  Wave  Height 

Figure  29  shows  the  daily  evolution  of  signifieant  wave  height  (H^)  in  the  SCS 
during  the  passage  of  TY  Muifa.  The  eolor  seale  symbolizes  the  value  of  ,  from  blue 

(0  m)  to  red  (16  m).  Here,  the  upper  panels  (Figure  29a-e)  are  for  the  period  before  TY 
Muifa  entering  the  SCS;  the  middle  panels  (Figure  29d-f)  are  for  the  period  during  TY 
Muifa  passing  through  the  SCS;  and  the  lower  panels  (Figure  29g-i)are  for  the  period 
after  TY  moving  out  from  the  SCS. 

Comparing  to  the  evolution  of  QTCWPM  wind  field  (Figure  11),  it  is  found  that 
the  distribution  of  follows  the  moving  TY  Muifa.  High  is  with  high  wind  speed 

and  high  waves  appear  to  the  right  of  the  typhoon  passage.  The  maximum  (16  m) 
oeeurs  as  TY  Muifa  reaehes  its  lowest  pressure  in  the  SCS  on  21  November. 

Another  high  oeeurs  in  east  of  Luzon  on  24  and  25  November  (Figure  29h-i). 

Comparison  between  upper  and  lower  panels  of  Figure  30  elearly  shows  the  differenee  in 
monsoon-generated  waves  before  TY  Muifa  arrival  (Figure  29a-e)  and  after  TY  Muifa 
departure  (Figure  29g-i).  This  indieates  that  the  existing  typhoon-generated  wave  field 
may  affect  the  generation  of  new  wind-wave  field. 

Figure  30a  shows  the  horizontal  distributions  of  WW3  generated  maximum 

and  maximum  wind  speed  during  TY  Muifa  passage  from  17  to  25  November.  The  core 
of  the  maximum  wind  speed  is  symmetric  along  the  track  of  the  TY  Muifa  passage 
(Figure  31b).  However,  the  core  of  the  maximum  is  asymmetric  in  width  and  in 

value  along  the  track  of  TY  Muifa  (Figure  30a)  with  higher  in  the  right  side  of  the 
track  than  the  left  side. 

TY  Muifa  has  its  lowest  pressure  at  Center-II  (11.6°A,  114. 4°F')  at  1800UTC  21 
November.  Figure  31  shows  the  horizontal  distribution  of  at  Center-II  with  the 
typhoon  translation  axis  being  rotated  towards  north.  The  radius  of  the  maximum  wind 
(7?„, )  is  14.0  km.  The  high  values  of  are  in  the  right  forward  quadrant  with  the 
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maximum  of  16  m  slightly  outside  of  from  the  typhoon  eenter.  This  eonsists  with 

earlier  studies  (Moon  et  al.  2003)  that  the  wave  generation  on  the  typhoon  eenter  has  an 
asymmetrie  distribution  in  four  quadrants. 

2.  Directional  Wave  Spectra  at  Center-II 

Figure  32  shows  the  daily  evolution  of  direetional  wave  speetra  at  Center-11.  The 
arrow  indieates  the  wind  speed  and  direetion  at  eaeh  moment.  The  winds  (7.8  m/s)  blow 
from  northeast  (45°  from  the  north)  before  1800UTC  20  November  and  after  1800UTC 
23  November,  showing  the  winter  monsoon  dominates.  In  these  two  periods,  a  set  of 
low-frequeney  waves  exists  and  aligns  with  the  wind  in  a  slight  angle  leftward. 
Conversely,  from  1800UTC  20  November  to  1800UTC  23  November,  rapid  ehange  of 
waves  in  both  direetion  and  frequeney  is  found.  Furthermore,  after  the  departure  of  TY 
Muifa,  evident  wave  paeket  still  remains  in  the  direetion  from  30  °  to  70  ° . 

To  investigate  variation  of  the  direetional  wave  speetra  during  TY  Muifa  passage, 
the  6-hourly  speetral  evolution  is  presented  (Figure  33)  from  1800UTC  20  November  to 
1800UTC  22  November.  It  is  seen  that,  while  the  wind  is  inereasing  and  turning 
eountereloekwise  as  TY  Muifa  passing  by,  the  waves  are  generated  along  with  the 
turning  wind  direetions.  As  the  typhoon  loeated  right  on  Center-11  at  1800UTC  21 
November,  the  wind  speed  signifieantly  deereases  and  the  wind  veetor  shifts  to  the 
opposite  direetion.  A  wide  range  of  waves  are  generated  with  monsoon  dominated  swell 
in  the  southwestward  direetion,  and  that,  as  typhoon  wind  ehanges  rapidly  in  speed  and 
direetion,  a  pair  of  wave  paekets  are  generated  at  the  time  of  typhoon  arrival  (Figure 
33e).  These  typhoon-generated  waves,  espeeially  in  the  opposite  direetion  to  monsoon 
winds  (i.e.,  the  northeastward  direetion)  deeay  fast  but  still  exist  four  days  after  the 
typhoon  departure. 

Four-way  loeations  are  seleeted  around  eaeh  eenter  with  the  maximum  wind 
radius  (/?„,)  away  to  investigate  the  effeet  of  translation  speed  on  the  direetional  wave 

speetra,:  (1)  forward  loeation  (F),  (2)  backward  location  (B),  (3)  rightward  location  (R), 
and  (4)  leftward  location  (L)  (see  up-left  panel  of  Figure  35  to  Figure  38).  The  typhoon 
center  is  indicated  by  ‘C’.  The  dashed  curves  are  circles  with  radius  from  the 
typhoon  centers.  The  hollow  arrow  indicates  the  typhoon’s  translation  velocity  vector. 
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The  solid  arrow  is  the  wind  vector  at  each  location.  Consistent  with  the  previous  study 
(Moon  et  al.  2003),  the  wave  spectra  in  the  forward  and  rightward  directions  show  that 
high-frequency  waves  are  generated  aligning  the  wind  direction  due  to  resonant  effects. 
But  the  wave  spectra  in  the  leftward  and  backward  directions  have  a  complicated 
structure. 

3.  Effect  of  Typhoon  Translation  Speed 

Figure  34  shows  the  6-hour  evolutions  of  at  the  typhoon  center  for  all  four 
points.  In  order  to  investigate  the  effect  of  typhoon  translation  speed,  all  the  fields 
have  bee  rotated  with  the  typhoon  translation  axis  pointing  toward  north.  Comparison  of 
the  fields  between  Center-I  and  Center-Ill  leads  to  the  effects  of  typhoon  translation 

on  the  wave  characteristics.  The  maximum  is  moving  faster  at  Center-I  than  at 
Center-Ill.  The  value  of  maximum  at  Center-I  increases  as  the  typhoon  low  pressure 
decreases.  It  is  found  that  the  highest  is  located  slightly  outside  away  from  the 
typhoon  center  in  the  right  forward  quadrant,  and  that,  the  higher  typhoon  translation 
speed  the  larger  angle  between  the  highest  and  the  eastward  direction  is. 

For  the  same  typhoon  intensity  (Center-I  and  Center-Ill),  the  directional  wave 
spectra  at  the  typhoon  center  (location-C)  is  more  asymmetrical  for  the  faster  translation 
speed  (Figure  35C)  than  for  the  slower  translation  speed  (Figure  37C).  High  wave 
energy  centers  occur  in  both  frontward  and  backward  directions  for  the  slower  translation 
speed. 

4.  Effect  of  Typhoon  Intensity 

Typhoon  has  the  same  translation  speed  but  different  intensity  at  Center-II  and 
Center-IV.  Comparison  of  Hs  and  directional  wave  spectra  between  Center-II  and 
Center-IV  leads  to  the  effect  of  typhoon  intensity.  Figure  34  shows  that  the  maximum 
is  at  least  5  m  larger  and  located  closer  to  the  typhoon  center  at  Center-II  than  at 

Center-IV.  This  is  because  the  typhoon  is  stronger  with  smaller  at  Center-II  than  at 
Center-IV.  As  the  typhoon  translation  speed  increases  significantly  at  Center-IV  (7.29 
m/s),  the  shape  of  contours  becomes  elliptic. 
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The  directional  wave  spectra  at  the  typhoon  center  clearly  show  the  following 
features.  The  waves  are  generated  in  both  frontward  and  backward  directions  for  higher 
typhoon  intensity  (Center-II,  Figure  36C),  but  only  in  the  frontward  direction  for  lower 
typhoon  intensity  (Center-IV,  Figure  38C).  Among  the  four  locations  around  the  typhoon 
center,  the  directional  wave  spectra  are  very  different  at  locations  (F,  R,  and  L)  but  less 
different  at  location  B  between  Center-II  and  Center-IV.  At  locations  (F,  R,  and  L),  two 
wave  packets  are  generated  in  both  upwind  and  downwind  directions  for  high  typhoon 
intensity  (Figure  36),  but  only  one  wave  packet  is  generated  in  downwind  direction  for 
low  typhoon  intensity  (Figure  38).  The  difference  between  Center-II  and  Center-IV 
(comparing  Figure  36  with  Figure  38)  is  much  larger  than  the  difference  between  Center- 
I  and  Center-Ill  (comparing  Figure  35  with  Figure  37).  This  indicates  that  the  effect  of 
typhoon  intensity  is  much  larger  than  the  effect  of  typhoon  translation  speed  on  the  SCS 
wave  characteristics. 

D.  EFFECTS  OF  MONSOON  WINDS 

When  the  ideal  typhoon  winds  during  TY  Muifa  passage  through  the  SCS  are 
used  as  the  wind  input  for  WW3,  the  effect  of  background  winds  (winter  monsoon)  is  not 
included.  Difference  between  two  simulations  described  in  this  and  previous  sections  is 
considered  the  effect  due  to  the  monsoon  forcing. 

1.  Significant  Wave  Height 

Figure  39  shows  the  daily  evolution  of  significant  wave  height  in  the  SCS 

during  the  passage  of  TY  Muifa.  Similar  to  Figure  29,  Figure  39a-c  are  for  the  period 
before  TY  Muifa  entering  the  SCS;  Figure  39d-f  are  for  the  period  during  TY  Muifa 
passing  through  the  SCS;  and  Figure  39g-i  are  for  the  period  after  TY  moving  out  from 
the  SCS. 

Comparing  Figure  39  to  Figure  29,  it  is  found  that  the  horizontal  distribution  of 
is  quite  similar  between  using  the  ideal  typhoon  winds  (Figure  39)  and  the 

QTCWPM  winds  (Figure  29)  near  the  typhoon  center.  But  the  maximum  using  the 

ideal  typhoon  winds  is  15  m  and  less  than  that  (16  m)  using  QTCWPM  winds.  Using  the 
ideal  typhoon  winds,  the  wave  energy  is  generated  near  the  typhoon  center  along  the 
typhoon  passage.  After  the  departures  of  TY  Muifa,  the  wave  energy  remains  for  one  to 
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two  days  and  then  disperses  along  the  typhoon  passage.  Wave  energy  propagates  to  the 
forward,  the  right-frontward,  and  the  left-rearward  direetions.  There  is  no  high  area 

oeeurring  east  of  Luzon  from  24  to  25  November. 

Figure  40  shows  the  horizontal  distributions  of  WW3  generated  maximum 

(Figure  40a)  and  maximum  wind  speed  (Figure  40b)  during  TY  Muifa  passage  using  the 
ideal  winds  from  17  to  25  November.  Similar  to  the  simulation  using  the  QTCWPM 
winds  (Figure  30),  the  eore  of  the  maximum  wind  speed  is  symmetrie  along  the  traek  of 
TY  Muifa  (Figure  40b).  However,  the  eore  of  the  maximum  is  asymmetrie  in  width 

and  in  value  along  the  traek  of  TY  Muifa  (eompared  to  Figure  40a)  with  higher  in  the 
right  side  of  the  traek  than  the  left  side. 

Comparing  Figure  40  with  Figure  30,  it  is  found  that  the  monsoon  wind  does  not 
affeet  Hs  near  the  typhoon  passage,  but  does  affeet  far  from  the  typhoon  traek, 
espeeially  in  the  northern  SCS.  In  the  maximum  wave  field,  the  maximum  field  is 

generally  higher  with  monsoon  present,  and  signifieant  difference  is  found  in  the  west  of 
Luzon. 

Figure  41  shows  mean  differences  of  and  wind  fields  between  the  two 

simulations.  The  difference  in  wind  fields  is  relatively  small  along  the  typhoon  passage 
and  relatively  large  away  from  the  typhoon  passage.  It  is  close  to  the  background  mean 
wind  (Figure  10).  The  maximum  difference  in  is  located  in  the  deep  water  area  away 

from  the  coast.  It  indicates  a  fetch-limited  case,  which  the  higher  waves  produced  farther 
from  the  coast  by  steady  monsoon.  The  obstruction  of  shoreline  and  small  islands  is  also 
seen  on  the  wave  propagation  along  the  monsoon  direction. 

2.  Directional  Wave  Spectra 

Figure  42  shows  the  daily  evolution  of  directional  wave  spectra  at  Center-II.  The 
arrow  indicates  the  wind  speed  and  direction  at  each  moment.  Distinction  between 
Figure  42  and  Figure  32  is  the  swell  generation  (no  swell  generation)  in  the 
southwestward  direction  with  QTCWPM  (ideal  typhoon)  winds.  This  difference  is  due  to 
the  background  monsoon  winds. 
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To  investigate  the  differenee  of  the  direetional  wave  speetra  between  using 
QTCWPM  and  ideal  typhoon  winds  during  TY  Muifa  passage,  a  6-hour  speetral 
evolution  is  presented  (Figure  43)  from  1800UTC  20  November  to  1800UTC  22 
November.  Comparing  Figure  43  to  Figure  33,  it  is  found  that  the  two  direetional  wave 
speetra  are  almost  the  same  during  0600UTC  21  November  to  1200UTC  22  November 
when  TY  Muifa  is  around  Center-II,  but  quite  different  for  other  time  periods.  Before 
TY  Muifa  arriving  Center-II  (I800UTC  21  November),  there  is  high  wave  energy  in  the 
southwestward  direetion  using  the  QTCWPM  winds  (Figure  33  a)  but  none  using  the  ideal 
typhoon  winds  (Figure  43a).  One  day  after  the  typhoon  moving  out  of  Center-II,  the 
typhoon-generated  waves  in  swell  direetion  dispersed,  and  the  swell  energy  appeared 
again  (see  Figure  33).  This  suggests  that  the  monsoon-generated  swell  does  not  deeay 
and  remain  through  the  typhoon  period. 

Figure  44  shows  the  distribution  and  direetional  wave  speetra  using  the  ideal 

typhoon  winds  on  the  typhoon  position  at  Center-IV  (water  depth:  116  m).  Figure  44  is 
almost  the  same  as  Figure  38,  whieh  indieates  that  even  when  the  typhoon  winds  weaken, 
the  typhoon  generated  waves  still  dominate. 

E.  EFFECTS  OF  BOTTOM  TOPOGRAPHY 

A  third  WW3  simulation  is  eonducted  for  the  SCS  with  everything  the  same  as 
Seetion  6  but  the  uniform  water  depth  (2000  m.  Figure  45).  The  surrounding  islands, 
Taiwan,  the  Philippians  Islands,  Palawan,  Borneo,  and  the  Indonesia  Islands,  have  been 
removed.  Differenee  between  the  third  and  the  seeond  simulations  leads  to  the  effeets  of 
island  obstruetion  and  bottom  topography  on  the  wave  eharaeteristies. 

1.  Significant  Wave  Height 

Figure  46  shows  the  daily  evolution  of  simulated  by  WW3  with  the  ideal 
typhoon  winds  with  uniform  water  depth.  The  effeets  of  topography  are  examined  by 
eomparing  Figure  46  to  Figure  40  (realistie  topography).  The  maximum  is  about 

15.5  m,  a  little  larger  than  using  the  realistie  topography  (15  m).  Similar  feature  is  also 
found  in  the  right-frontward  direetion  of  the  typhoon  eenter,  expeet  the  right-front 
quadrant  has  higher  value  and  wider  distribution.  As  TY  Muifa  enters  the  SCS  on  20 
November,  the  typhoon-generated  waves  are  well-developed  due  to  no  island  obstruetion. 
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The  wave  propagations  in  the  right-rearward  direction  of  the  typhoon  passage  are  similar, 
but  in  the  uniform  water  depth  simulation,  it  is  much  clear  that  the  wave  does  propagate 
forward,  and  faster  than  the  typhoon  translation.  Although  there  is  no  island  blockage  in 
the  left-rearward  direction,  the  wave  propagations  are  still  limited. 

It  is  evident  that  along  the  typhoon  passage,  the  typhoon  generated  waves  more 
like  to  propagate  in  the  right-backward  and  left-forward  directions,  and  the  right-forward 
and  left-backward  directions  are  restricted.  This  may  be  due  to  the  momentum  flux  into 
the  typhoon  passage  in  the  right-forward  and  left-backward  quadrants. 

Figure  47  shows  the  maximum  values  of  the  simulated  waves,  and  the  mean 
difference  between  normal  water  depth  and  the  uniform  water  depth.  Compared  Figure 
47a  with  Figure  40a,  it  is  found  that  higher  appears  from  \  \9° E  to  \22° E ,  and  that, 

as  TY  Muifa  close  to  the  Asian  landmass,  the  wave  energy  does  not  disperse  due  to  no 
change  in  water  depth,  and  disappear  only  when  the  waves  contact  coastline.  The  depth 
of  the  nearshore  dispersion  area  is  between  20  to  50  m. 

This  feature  is  exhibited  better  in  the  mean  7/^  difference  between  normal  depth 
and  uniform  depth  (Figure  47b).  The  color  scale  is  from  0  m  (blue)  to  1.5  m  (red).  It 
shows  that  from  the  uniform  depth  is  developed  better  than  from  the  real  SCS 
topography,  and  the  maximum  difference  is  within  1.5  m.  Large  difference  is  seen  along 
the  typhoon  passage  while  close  to  coastline.  The  notable  difference  is  also  seen  in  the 
Gulf  of  Thailand,  around  Palawan,  and  in  the  west  water  of  Borneo.  In  the  east  water  of 
Borneo,  previous  simulation  suggests  that  no  wave  propagation  into  it,  but  in  this  uniform 
water  simulation,  the  wave  propagation  is  observed. 

It  is  observed  that  the  SCS  wave  does  not  propagate  into  the  Sulu  Sea  due  to  the 
obstruction  of  Palawan  and  Borneo.  Therefore,  the  Sulu  Sea  may  be  considered  as  a 
different  wave  system  from  the  SCS,  and  the  wave  simulation  of  the  Sulu  Sea  in  this 
study  may  be  underestimated.  Compared  the  difference  in  the  south-west  SCS  with  the 
bathymetry  contour  of  the  SCS  (Figure  24),the  difference  of  0.5  m  does  follow  100  m 
water  depth  properly. 
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2.  Directional  Wave  Spectra 

Figure  48  shows  the  daily  evolution  of  direetional  wave  speetra  with  the  idealized 
typhoon  winds  and  uniform  water  depth  on  the  same  loeation  as  Figure  32.  Compared 
Figure  48  with  Figure  32,  it  is  seen  that  few  wave  energy  reaehes  before  the  typhoon 
arrival,  but  after  the  typhoon  departure,  the  wave  speetra  are  identieal.  It  may  due  to  the 
wave  energy  does  not  disperse  by  the  Philippians  islands  and  propagate  here  before  TY 
Muifa  arrived. 

Figure  49  shows  the  distribution  and  direetional  wave  speetra  using  the  ideal 

typhoon  winds  and  uniform  water  depth  on  the  typhoon  position  at  Center-IV. 
Comparing  Figure  49  (uniform  water  depth)  to  Figure  44  (realistie  topography,  water 
depth:  116  m),  it  is  found  that  the  direetional  wave  speetra  and  are  almost  the  same, 

whieh  indieates  that  the  topographie  effeet  on  the  wave  eharaeteristies  is  not  important 
when  the  loeal  typhoon  winds  dominates  with  the  water  depth  deeper  than  100  m. 
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Figure  24  The  SCS  bathymetry  from  the  TerrainBase  dataset.  Red  line  and  points  indicate 

the  passage  of  TY  Muifa  (2004). 
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Figure  25  TOPEX/Poseidon  crossover  points  in  the  SCS. 


Figure  26  WW3  result  compared  with  TOPEX/Poseidon  observations  on  all  crossover 

points  during  the  period  of  TY  Muifa  (2004).  (a)  The  paired  data  distribution,  and 
(b)  the  histogram  distribution  of  difference. 
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Figure  27  (a)  The  translation  speed,  (b)  the  lowest  pressure,  and  (c)  the  maximum  wind 

speed  of  TY  Muifa  (2004)  with  (d)  the  water  depth  along  the  typhoon  passage. 
Red  lines  indicate  four  designed  typhoon  centers. 
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Figure  29  Daily  evolution  of  in  the  SCS  during  the  period  of  TY  Muifa  (2004). 
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Figure  30  Total  distribution  of  (a)  the  maximum  from  WW3,  and  (b)  the  maximum 

QTCWPM  wind  speed  during  the  passage  of  TY  Muifa  (2004). 
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FigureSl  of  the  lowest  pressure  center  of  TY  Muifa  in  the  SCS  on  (1 1.6°A^ ,  \\AA°E) 

at  1800UT  21  November.  The  estimated  maximum  wind  radius  is  about  14.0  km. 
The  translation  direction  (hollow  arrow)  is  rotated  to  northward. 
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Figure  32  Daily  evolution  of  directional  wave  spectra  with  QTCWPM  winds  on  designed 
Center-II  (11.6°A^  and  114. 4°£’).  TY  Muifa  arrived  this  point  at  1800UTC  21 
November.  The  arrow  presents  QTCWPM  wind  speed  and  direction.  The  wind 
speed  value  is  100  times  of  the  axis  scale. 
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Figure  33  Detail  evolution  of  Figure  32  with  time  step  6  hours  from  1800UTC  20  November 

to  1800UTC  21  November. 
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Figure  34  6-hour  evolutions  of  on  four  designed  typhoon  centers  (a)  before  typhoon 

arrivals  (left  column),  (b)  typhoon  right  on  centers  (central),  and  (c)  after  typhoon 
departures  (right).  The  translation  directions  are  rotated  to  northward. 

72 


Frequency  (Hz)  Frequency  (Hz)  Distance  from  TY  Center  (km) 


100 


SWH  11/21  00:00:00 


(C)enter  of  TY 


(F)otward 


Frequency  (Hz) 


(R)ightwarcl  R^^ 


(L)ellward  R^ 


(B)acl<ward  R^^^ 


Frequency  (Hz) 


Frequency  (Hz) 


Figure  35  H ^  of  Center-I  (OOOOUTC  21  November,  11.9°A^  and  117.2°£'),  and  its  eentral 

(C)  and  four-way  (F,  R,  L,  and  B)  directional  wave  spectra.  The  hollow  arrow  is 
the  translation  direction.  The  format  of  wave  spectra  is  similar  to  Figure  32. 
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Figure  36  ,  and  central  (C)  and  four-way  (F,  R,  L,  and  B)  directional  wave  spectra  of 

Center-II  (1800UTC  21  November,  11.6°A^  and  114.4°£').  Similar  format  as 

Figure  35. 
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Figure  37  ,  and  central  (C)  and  four-way  (F,  R,  L,  and  B)  directional  wave  spectra  of 

Center-Ill  (1800UTC  22  November,  10.5  ° N  and  1 12.1  ° E).  Similar  format  as 

Figure  35. 
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Figure  38  ,  and  central  (C)  and  four-way  (F,  R,  L,  and  B)  directional  wave  spectra  of 

Center-IV  (0600UTC  24  November,  8.8  and  108.8  °£').  Similar  format  as 

Figure  35. 
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Figure  39  Daily  evolution  of  from  ideal  typhoon  winds  in  the  SCS  during  the  period  of 

TY  Muifa  (2004). 
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Figure  40  Total  distribution  of  (a)  the  maximum  and  (b)  the  maximum  wind  speed  from 

ideal  typhoon  winds.  Similar  format  as  Figure  30. 


Figure  41  (a)  The  mean  difference  of  between  Figure  29  and  Figure  39,  and  (b)  the 

mean  difference  of  wind  speed  between  Figure  1 1  and  Figure  13. 
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Figure  42  Similar  format  as  Figure  32  on  the  same  location  (1 1.6°  and  1 14.4  “£’)  with 
ideal  typhoon  winds.  The  time  step  is  24  hours  from  1800UTC  16  November  to 

1800UTC  24  November. 
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Figure  43  Detail  evolution  of  Figure  42  with  time  step  6  hours  from  1800UTC  20  November 

to  1800UTC  21  November. 


80 


Frequency  (Hz)  Frequency  (Hz)  Distance  from  TY  Center  (km) 


SWH  11/24  06:00:00 


(F)otward 


Frequency  (Hz) 


Frequency  (Hz) 


Frequency  (Hz) 

(B)ackward 


Figure  44  ,  and  central  (C)  and  four-way  (F,  R,  L,  and  B)  directional  wave  spectra  from 

ideal  typhoon  winds  of  Center-IV  (0600UTC  24  November,  8.8  °A^  and 
108.8  ° E).  Similar  format  as  Figure  35. 
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Figure  46  Daily  evolution  of  from  ideal  typhoon  winds  with  uniform  depth  in  the  SCS 

during  the  period  of  TY  Muifa  (2004). 
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Figure  47  (a)  The  total  distribution  of  the  maximum  from  ideal  typhoon  winds  with 

uniform  water  depth.  Similar  format  as  Figure  40a.  (b)  The  mean  difference  of 
between  Figure  39  and  Figure  46.  Similar  format  as  Figure  41a. 
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Figure  48  Similar  format  as  Figure  32  on  the  same  location  (1 1.6°  and  1 14.4  “£’)  with 
ideal  typhoon  winds  and  uniform  water  depth.  The  time  step  is  24  hours  from 
1800UTC  16  November  to  1800UTC  24  November. 
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Figure  49  ,  and  central  (C)  and  four-way  (F,  R,  L,  and  B)  directional  wave  spectra  of 

Center-IV  (0600UTC  24  November,  S.S° N  and  108.8 °E)  from  ideal  typhoon 
winds  with  uniform  water  depth.  Similar  format  as  Figure  35. 
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Switeh  Parameters 

Charaeteristies 

DUM 

Dummy  to  be  used  if  WW3  is  to  be  installed  on  previously 
untried  hardware 

LRB8 

8  byte  words 

SHRD 

Shared  memory  model,  no  message  passing 

SEED 

Seeding  of  high-frequeney  energy 

GRDl 

Settings  direetly  hardwired  to  user-defined  spatial  grids 
(spherieal  eoordinate  with  1/4°  grids) 

SPl 

Eiser-defmed  speetral  grids. 

PR2 

QETCKEST  EiETIMATE  propagation  seheme  with  Booij  and 
Holthuijsen  (1987)  dispersion  eorreetion 

ST2 

Tolman  and  Chalikov  (1996)  souree  term  paekage 

STAB2 

Enable  stability  eorreetion  for  Tolman  and  Chalikov  (1996) 
souree  term  paekage 

NEl 

Nonlinear  interaetion  (DIA) 

BTl 

JONS  WAP  bottom  frietion  formulation 

WIND2 

Approximately  quadratie  interpolation 

CUR2 

Approximately  quadratie  interpolation 

ol 

Output  of  boundary  points  in  grid  preproeessor 

o2 

Output  of  the  grid  point  status  map  in  grid  preproeessor 

o2a 

Generation  of  land-sea  mask  file  mask.WW3  in  grid 
preproeessor 

o3 

Additional  output  in  loop  over  fields  in  field  preproeessor 

o4 

Print  plot  of  normalized  1-D  energy  speetrum  in  initial 
eonditions  program 

o5 

2-D  energy  speetrum 

o6 

Spatial  distribution  of  wave  heights  (not  adapted  for  distributed 
memory) 

o7 

Eeho  input  data  for  homogeneous  fields  in  generie  shell 

Table  4  The  WW3  model  switeh  parameters  setting. 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  investigated  the  SCS  wave  charaeteristies  using  WW3  with  winter 
typhoon  forcing.  The  model  was  forced  by  high  resolution  wind  field,  which  was 
obtained  from  the  QuikSCAT  observation  combined  with  the  Tropical  Cyclone  Wind 
Profile  Model  (TCWPM).  The  WW3  model  was  evaluated  using  the  TOPEX/Poseidon 
altimetry  observation  during  the  period  of  Typhoon  Muifa  (2004). 

The  effects  of  typhoon  and  monsoon  on  the  wave  characteristics  were  analyzed 
with  three  scenarios  using;  (1)  realistic  winds  and  bottom  topography,  (2)  ideal  typhoon 
winds  and  realistic  bottom  topography,  and  (3)  ideal  typhoon  winds  and  uniform  bottom 
topography.  The  typhoon  effect  was  simulated  using  the  ideal  typhoon  wind  (Scenario- 
2).  The  monsoon  effect  was  simulated  using  the  difference  between  Scenario- 1  and 
Scenario-2.  Use  of  the  uniform  bottom  topography  represented  the  open  ocean  case.  The 
topographic  effect  was  examined  by  the  difference  between  Scenario-2  and  Scnario-3. 

Along  the  typhoon  translation  track,  the  core  of  the  maximum  significant  wave 
height  (H^)  was  asymmetric  with  higher  and  wider  core  in  the  right  side  than  the  left 
side.  The  maximum  was  about  16  m  located  at  minimum  typhoon  lowest  pressure 
center.  At  a  single  position,  the  maximum  was  always  in  the  right-front  quadrant  of 
the  typhoon  center  and  located  slightly  outside  the  maximum  wind  radius  ( ).  Before 

typhoon  arrival,  the  monsoon-generated  swell  dominated  and  maximum  significant  wave 
height  was  about  4  m;  after  typhoon  departure,  the  monsoon  swell  remained  dominant 
and  the  typhoon-introduced  waves  decayed  gradually  and  lasted  for  three  days. 

The  waves  were  generated  in  both  frontward  and  backward  directions  for  higher 
typhoon  intensity,  but  only  in  the  frontward  direction  for  lower  typhoon  intensity.  When 
the  typhoon  intensity  was  reduced  (from  Center-II  to  Center-IV),  the  directional  wave 
spectra  were  very  different  at  the  frontward,  rightward,  and  leftward  locations.  Two 
wave  packets  were  generated  in  both  upwind  and  downwind  directions  for  high  typhoon 
intensity,  but  only  one  wave  packet  was  generated  in  downwind  direction  for  low 
typhoon  intensity.  However,  the  directional  wave  spectra  were  less  different  at  the 
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backward  location.  For  the  typhoon-generated  waves  in  the  SCS,  the  faetors  affeeting  the 
SCS  wave  eharaeteristies  with  their  importanee  from  large  to  small  are  listed  as  follows: 
typhoon  intensity,  distanee  and  direetion  from  typhoon  eenter,  the  typhoon  translation 
speed,  monsoon,  and  topography. 

The  typhoon  intensity  had  a  greater  effeet  on  the  SCS  waves  than  the  typhoon 
translation  speed.  As  the  typhoon  passed  by,  wind  speed  and  direetion  ehanged  rapidly. 
Stronger  intensity  with  slower  translation  speed  usually  generated  a  pair  of  wave  paekets. 
On  the  eontrary,  typhoon  foreing  had  greater  effeet  on  the  SCS  waves  than  monsoon 
foreing.  Low-frequeney  swell  was  generated  by  monsoon  foreing.  When  the  typhoon 
was  present,  the  typhoon-introdueed  waves  dominated  and  were  mueh  greater  than 
monsoon-generated  swell. 

In  wave  propagation,  the  typhoon-generated  wave  transformed  aeross  the  SCS  to 
farther  than  1000  km  but  only  in  the  forward,  the  left-frontward,  and  the  right-rearward 
direetions  of  the  typhoon  eenter.  The  wave  energy  did  not  propagate  to  the  left-rearward 
direetion  of  the  typhoon  eenter.  Along  the  wave  propagation  direetion,  the  obstruetions 
of  eoastline  and  small  islands  generated  the  wave  shadow  zones.  This  feature  implies 
that  the  SCS,  a  semi-elose  marginal  sea,  ean  be  eonsidered  as  an  independent  wave 
system  from  nearby  oeeans  exeept  the  energy  transform  through  the  Luzon  Strait,  whieh 
is  deep  to  1800  m.  The  Sulu  Sea,  whieh  separated  by  the  Palawan  Island,  ean  be  seen  as 
a  different  wave  system  from  the  SCS. 

The  bathymetrie  effeet  was  important  only  for  water  depth  shallower  than  100  m. 
For  depth  deeper  than  100  m,  the  bathymetry  may  have  an  effeet  on  ,  but  almost  no 
effeet  on  the  wave  propagation  direetion  and  deeay  frequeney. 

The  field  observation  in  tropieal  eyelones  of  Powell  et  al.  (2003)  reports  that  the 
drag  eoeffieient  ( )  is  mueh  less  than  previously  thought  in  extreme  high  wind  speed, 

i.e.,  above  40  m/s.  At  high  wind  speed,  the  overestimate  of  the  drag  eoeffieient  (C^)  in 

WW3  has  also  been  reported  (Moon  et  al.  2004a).  One  possible  improvement  is  to  use  a 
new  seheme  to  eompute  the  drag  eoeffieient.  Furthermore,  this  high  wind  speed 
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overestimate  may  be  examined  by  comparing  wave  simulation  with  field  wave 
observation,  if  detail  wave  measurements  (e.g.,  motion  buoys)  are  available  during  the 
typhoon  passage  in  the  future. 
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